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LONGITUDINAL  STABILITY  AND  CONTROL 
OF  A  DUCTED-ROTOR  FLYING  PLATFORM 


SUMMARY 


An  analysis  is  presented  for  the  prediction  of  the  longitudinal 
stability  and  control  characteristics  of  a  ducted -rotor .Flying  Platform. 
The  stability  characteristics  of  this  vehicle  are  corroborated  by  the 
investigations  and  exhibited  rather  clearly  by  various  stability  diagrams. 
Artificial  stabilization  is  discussed  and  can  be  used  to  stabilize 
certain  configurations.  The  predictions  of  excessive  trim  requirements 
and  large  tilt  angles  for  forward  flight  are  substantiated.  Also,  the 
effects  of  changing  duct  design  and  vertical  center  of  gravity  location 
on  the  stability  derivatives  are  discussed,  and  material  is  presented 
to  reflect  the  degree  to  which  these  changes  can  be  used  as  aids 
in  stabilizing  the  Platform. _ 
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LONGITUDINAL  STABILITY  AND  CONTROL 
OF  A  DUCTED -ROTOR  FLYING  PLATFORM 

INTRODUCTION 

The  ducted-rotor  Flying  Platform  is  essentially  a  stand-on  type,  VTOL, 
lifting  device  operated  by  a  pilot.  Lift  is  produced  by  two,  counter-rotating, 
non-articulated ,  rotors  located  below  the  pilot's  stand-on  platform  and 
powered  by  reciprocating  engines.  (Fig.  1) 

The  conventional  control  systems  of  airplanes  and  helicopters  are 
eliminated,  since  control  is  primarily  obtained  by  the  method  often  refer¬ 
red  to  as  "kinesthetic  control".  (Ref.  1) 

Longitudinal  and  lateral  control  is  maintained  when  the  pilot,  acting 
as  an  automatic  feedback  by  his  instinctive  balancing  reactions,  provides 
the  correcting  moment  necessary  to  control  and  stabilize  the  machine  after 
a  disturbance.  In  steady  forward  flight  he  provides  the  required  control 

moment  by  shifting  his  weight  and  leaning  in  the  direction  of  the  desired 

£ 

motion. 

Directional  control  is  obtained  through  the  movement  of  yaw  vanes 
located  in  the  duct  and  operated  by  the  twist  grip  handle  bars,  which  also 
provide  a  method  for  power  adjustment. 

In  this  Report,  all  reference?  to  the  term  '  Flying  Plationn"  snould 
be  construed  to  mean  the  subject  vehicle. 

The  purpose  of  this  Report  is  fo  predict  and  understand,  by 
theoretical  and  ..^lytical  studies,  the  attainable  flying  qualities  asso¬ 
ciated  with  the  stability  and  control  of  the  ducted-rotor  Flying  Platform 
as  affected  by  changes  in  duct  design,  vertical  canter  of  gravity  location 
and  artificial  stabilization. 
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The  ducted-rotor  Flying  Platform  has  been  flying  for  many  years 
and  several  models  have  been  built. 

Stability  and  control  problems  have  developed  for  the  longitudinal 
mode  of  motion  of  these  machines.  The  main  difficulties  are  associated 
with  large  nose-up  pitching  moments,  poor  gust  response,  large  tilt 
angles  for  forward  flight  and  marginal  handling  qualities . 

For  the  studies  in  this  Report,  equations  of  motion  based  on  small 
perturbation  theory  were  developed  and  the  static  stability  derivatives 
were  calculated  from  data  obtained  by  wind  tunnel  tests  performed  at 
the  David  Taylor  Model  Basin.  (Ref.  2) 

Variations  in  vertical  center  of  gravity  location,  together  with  radi 
cally  different  duct  designs,  (Fig.  2),  were  analyzed  for  their  effect  on 
the  stability  and  control  of  the  vehicle. 

The  effects  of  artificial  stabilization  of  the  type  used  on  vehicles 
of  this  design  were  also  considered. 
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NOTATION 
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A,  B ,  C ,  D ,  E 

A  .  B  c  ,  D 
TI  H'  H'  H 


slope  of  blade  section  lift  curve  per  radian 
number  of  blades 

o 

chord  length  of  duct  (ft.) 
equivalent  flat  plate  area  of  pilot 
(approximately  9  sq.  ft.) 

non-dimensional  moment  of  inertia  parameter  (0.0108) 
radius  of  gyration  (ft.) 

distance  between  drag  center  of  pilot  and  center 
of  gravity  of  gross  weight  aircraft  (ft.) 
mass  of  complete  aircraft  (slugs)  ' 
poles  for  root  locus  for  C 

m 

dynamic  pressure  (lbs.  per  sq.  ft.) 
time  (sec.) 

distance  between  the  of  action  of  thrust  force  and 
center  of  gravity  of  gross  weight  aircraft  (ft.) 
zero  for  root  locus  for  C 

m 

distance  between  line  of  action  of  pilot  drag  force 
and  center  of  gravity  of  gross  weight  aircraft 
(ft.) 

distance  between  line  of  action  of  H-force  and 
center  of  gravity  of  gross  weight  aircraft  (ft.) 
coefficients  of  the  characteristic  equation  in 
forward  flight 

coefficients  of  the  characteristic  equation  in 
hover 
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Ci,CrC3,DrD2/D3  increments  to  coefficients  of  characteristic 


£1'E2 


,  2 
*  j  ^ 


equations  for  C  and  C  in  root  locus  equation 
mm 

M  * 


CD  =  CD>/2 


drag  coefficient  of  pilot  where  c  •  =  f/S 

D 


H 


'R 


H  pTTR2  (  flR)2 


H-force  coefficient 


m  ’  prrR^  ft  R)2r 


pitching  moment  coefficient 


Cm  = 


T  prrR2  (£2R)2 


thrust  coefficient 


'W  p  -rr  R^ 


¥■ 


weight  coefficient 


acT 


O' 


d£j_ 


dn 


H 


dpH 


dcH 


'H 


dM 


m 


a 


6a 

3  Cm. 


m 


dM 


c 

c 


_  3ch 

Hd(Ao)  W3J 

=  _ 

md (  A 9)  ^d(ie» 


static  stability  force  and  moment 
coefficient  derivatives ,  (moments 
taken  about  the  quarter  chord) 


CT  '  °T  '  CH  '  °H  ' 
ot  n  a  a 


c  ,  c  ,  c 

m  m  m  ,  .  a  . 
a  p  d  ( A  6) 


overall  values  of  force  and  moment  coefficient 
derivatives 


D  =  fp  V  /  2 
F. 


drag  force  of  the  pilot  (lbs.) 


‘X 


horizontal  force  (lbs) 

total  force  measured  along  X  body  axis  (lbs.) 
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H 

Iy  =  mky2 


It 

L.C.D. 

9 

M  =  M,,  -  M 


M 


c/4 


mp 

My 

eg 

P 


R 

$LBP 

V"r! 

S.C.D.. 

T 


V 

w 


total  force  measured  along  Z  body  axis  (lbs.) 

total  force  parallel  to  rotor  disc  (lbs.) 

pitching  moment  o t  inertia  of  gross  weight 

2. 

aircraft  (approximately  125  slug-ft.  ) 
gain  of  artificial  stabilization  feedback  loop 
Long  Chord  Duct 

pitching  moment  excluding  pilot  control  moment, 
measured  about  Y  axis  (ft. -lbs.) 
aerodynamic  moment  created  by  the  aircraft, 
without  pilot,  measured  about  the  c/4 
point  (ft. -lbs.) 

control  moment  exerted  by  the  pilot  (ft. -lbs.) 

total  moment  measured  about  the  Y  axis  (ft. -lbs.) 

* 

period  of  oscillation  (sec  ) 

also  poles  for  root  locus  for  C 

a 

rotor  radius  (3 . 5  ft. ) 

also  Routh's  discriminant 

abbreviation  for  Root  Locus  Break  Points 
2 

area  of  rotor  (ft.  ) 

e 

Short  Chord  Duct 

total  thrust,  force  of  duct  and  rotor  measured 
normal  to  the  rotor  disc  (lbs.) 

*lso  time  to  half  or  double  amplitude  (sec.) 
also  time  constant  for  phase  lag  of  artificial 
stabilization  feedback  loop 
resultant  velocity  of  aircraft  (ft.  per  sec.) 
gross  weight  of  complete  aircraft  including 
pilot  (550  lbs.) 


X,Y,Z 

Z 

a 


7 

9 


P 

m 

T=  piTR^TR 
O 


General 


A(  ) 


6 

body  a$is  system 
zero  for  root  locus  for  C 

m 

a 

angle  of  attack,  angle  between  relative^ind 
vector  and  plane  perpendicular  to  the  rotor 
axis  (radians) 

flight  path  angle  measured  with  reference  to 
horizon  (radians) 

aircraft  pitch  angle  with  reference  to  horizon 
(radians) 

root  of  the  characteristic  equation 
tip  speed  ratio  or  velocity  ratio 

air  density  (slugs  per  cubic  ft.) 
aerodynamic  time  £sec.) 

rotor  angular  velocity  (radians  per  sec.) 

*  »  * 

differential  operator  forms  for  time  ratio  t/r 

differential  operator  forms  for  time  t 

subscript  zero  designates  initial  condition 
before  perturbation 
designates  perturbation  quantities 
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DEVELOPMENT 

Equations  of  Longitudinal  Motion 

The  development  of  the  equations  of  motion  for  this  system  Is  based 
on  the  commonly  used  methods  for  the  study  of  aeronautical  dynamics. 

The  basic  theory  is  the  familiar  small  perturbation  theory  that  assumes 
the  disturbed  motion  of  the  aircraft  to  be  one  of  small  oscillations  about 
some  steady-state  flight  condition. 

As  usual,  the  equations  of  motion  are  written  in  accordance  with 
the  Newtonian  laws  of  motion  with  the  important  requirement  that  all 
accelerations  and  rates  of  change  of  moments  of  momentum  must  be 
expressed  along  alkes  fixed  in  space. 

Equations  of  motion,  as  well  as  other  important  work  in  this 
area,  have  been  developed  at  the  Advanced  Research  Division  of  Hiller 
Aircraft  Company  by  G.  J.  S^ssingh  and  A.  H.  Sacks. 

Also,  dynamic  stability  of  a  machine  similar  to  a  Flying  Platform 
has  been  treated  at  Princeton  University  in  Refs.  3  and  4. 

The  development  of  the  equations  of  motion  in  this  Report  follows 
the  latter  work,  the  principle  differences  being  in  definition  and  notation. 

For  the  specific  case  of  the  longitudinal  motion  of  the  Flying 
Platform,  we  have  a  dynamic  system  with  three  degrees  of  freedom;  namely, 
translation  forward  or  backward,  translation  upward  or  downward  and 
rotation  about  the  pitch  axis . 

This  system  will  be  described  sufficiently  by  three  equations  of 
motion  for  the  so-called  control-fixed  case  similar  to  that  studied  in 
airplane  dynamics.  The  control-fixed  case  of  the  Flying  Platform  can 
be  thought  of  as  the  condition  where  the  pilot  remains  rigid  on  the 
Platform,  from  the  initial  condition,  and  induces  no  new  or  additional 
increments  of  control  movement  throughout  the  motion  to  be  studied. 


ft******* 
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In  order  to  facilitate  the  description  of  these  motions,  choice  is 
made  of  the  Eulerian  or  moving  axis  system  fixed  within  the  aircraft  and 
always  moving  with  it.  For  this  axis  system,  the  standard  N.A.S.A. 
right-hand  axis  system  is  selected  with  the  positive  X,  Y,  Z,  axes 
in  the  forwaro,  right  side,  and  downward  directions  respectively. 
Choosing  these  axes  as  body  axes ,  the  origin  will  always  be  located  at 
the  center  of  gravity  and  will  move  slightly  with  respect  to  the  aircraft 
whenever  the  center  of  gravity  is  shifted.  The  positive  X  axis  is  directed 
toward  the  front  of  the  aircraft,  parallel  to  the  rotor  plane,  and 
aligned  with  the  relative  wind  in  azimuth  but  not  necessarily  in  pitch. 

The  positive  Z  axis  is  directed  downward,  perpendicular  to  the  rotor 
plane,  i.e.  aligned  parallel  to  the  thrust  line  or  rotor  shaft. 

The  Y  axis  will  be  the  aircraft  pitching  axis  located  at  the 
center  of  gravity.  This  axis  system,  with  notation  for  the  forces, 
moments,  and  displacements,  is  shown  in  Fig.  3. 

In  this  discussion,  the  forces  and  moments  to  be  resolved  along 
and  about  the  body  axis  system,  together  with  other  notation,  are 
described  as  follows: 

My  total  moment  measured  about  the  Y  axis  located  at  the 

c .  g . 

center  of  gravity.  Clockwise  or  nose-up  moments 
are  considered  positive. 

Mc/4  aerodynamic  moment  created  by  the  aircraft  without  pilot 

and  measured  positive  clockwis&-~(nose-up)  about 
the  quarter  chord  point.  (This  notation  arises  from 
the  wind  tunnel  balance  design  for  measuring  H- 
force  and  moments  in  connection  with  data  re¬ 
duction  when  transferring  moments  to  an  arbitrary 


axis  located  at  the  so-called  quarter  chord  point. 

This  quarter  chord  axis  is  parallel  to  the  Y  axis 
and  passes  through  the  thrust  line  at  the  quarter 
chord  point.) 

total  gross  weight  of  the  aircraft,  including  pilot  and 
all  components  „ 

force  measured  positive  upward  along  rotor  shaft 

line  .  perpendicular  to  the  rotor  plane  and  parallel, 
but  directed  opposite  to  the  positive  Z  direction, 
force  measured  perpendicularly  to  T  force  and  in  the  X 
direction.  It  is  taken  as  acting  in  the  quarter 
chord  plane  parallel  to  the  rotor  plane.  The  H- 
force  is  positive  (as  in  the  case  of  forward  flight) 
when  acting  aft  or  in  the  negative  X  direction, 
drag  of  the  pilot  measured  parallel  to  the  relative  wind 
and  acting  at  the  d/  ig  center  of  the  pilot.  Drag 
Center  and  center  of  gravity  location  of  the  pilot 
alone  were  assumed  to  coincide.  Also,  aerodynamic 
lift  caused  by  the  pilot  was  assumed  to  be  neg¬ 
ligible  . 

moment  arm  for  T  force  or  displacement  of  center  of  gravity 
from  thrust  line  measured  along  X  direction.  Dis¬ 
tance  is  positive  when  center  of  gravity  is  displaced 
aft  of  thrust  line,  and  negative  when  forward  of 
thrust  line  as  in  the  case  of  forward  flight, 
moment  arm  for  drag  force  of  pilot  where  1  cosa 

and  lp  is  the  distance  from  the  total  center  of 
gravity  position  to  the  drag  center  of  the  pilot 


as  shown  in  Fig .  1  . 


m 


2  moment  arm  for  H-force  or  distance  along  Z 

ri 

direction,  between  the  center  of  gravity  and 
the  line  of  H-force. 

a  angle  of  attack,  measured  between  the  rotor  tip  path 

plane  and  the  relative  wind  where  nose-up  is 
positive  alpha.  (Alpha  is  negative  in  forward  flight.) 

0  attitude  or  pitch  angle,  measured  between  rotor  tip  path 

plane  and  the  horizon  where  nose-up  is  positive 
tl  eta. 

7  flight  path  angle  measured  to  the  horizon.  As  shown  in 

Fin.  3,  7  =  0  -  a  where  gamma  is  negative  as  shown. 

In  the  development  Df  the  equations  of  motion  of  the  Flying  Platform, 
it  is  assumed  that  the  disturbed  motion  of  the  aircraft  is  one  of  small 
displacements  or  oscillations  about  each  steady-state  equilibrium-trim 
flight  condition.  This  steady-state  equilibrium-trim  condition  is  achieved 
whenever  the  summation  of  all  the  aerodynamic  and  gravity  forces  acting 
on  the  aircraft  equals  zero.  Also,  with  regard  to  the  inertia  terms,  the 
Newtonian  laws  of  motion  apply  to  any  axis  system,  provided  that  the 
forces  and  moments  are  resolved  along  and  about  the  particular  axis 
system  and  that  the  accelerations  expressed  along  and  about  these 
axes  are  the  true  accelerations  with  reference  to  fixed  space. 

With  these  considerations  the  dynamic  equations  of  motion  may  be 
written  in  the  form: 


IK-  E 


AFROO  YNA  M 1C 


EF  =  F 


^AERCL  YNAMIC 


EM.  =|1 


Xr  n 


Hf- 

Ty 

+  fy  =° 

(la) 

'  6#  AV/ry 

'  INERTIA 

f 

L 

+  r  -  o 

(lb) 

*  Gravity 

2  INERTIA 

-h 

M..  =o 

(lc) 

'aerodynamic 


c.g. 


INERTIA 


il 

Resolving  the  aerodynamic  and  gravity  contributions  along  body 
axes  (Fig.  3),  the  steady-state  equations  of.  motion  become: 

E  +  E  -  ~  Dccso£-  sin  6  =  0  (2a) 

XA 

Fh  ^  =  -T-  D  Sin  oC  COS  0  ~  0  (2b) 

A  G 

f  K  =■  fl  c/+-  H  •  hh + T*  xT+ D  ■  hd=  o  <20 

**  A 

Using  the  small  angle  assumptions  where: 

Sift  of  ~  c/"  COSof  «  / 

(3) 

sin  6  ~  6  cos  6  =5?  / 

the  steady-state  equations  are: 

Fx  +•  F  =  -  H  -  [)  -  W  Q  =  o  wa) 

p£  t  F  =  0  (4b) 

A  2  $ 

MYc.g*=  Mc/+-  H'Zh  +T'Xt  +  D'*d  "0  (40 

Expressing  the  inertia  contributions  by  using  the  Newtonian 
frame  of  reference  of  fixed  axes  and  body  axes  as  the  moving  axes  , 
the  inertia  terms  become: 

#  •  #  • 

F,  =  - /m.  V  cos  of  ~/rn  V  Y  siK  c/T  =  -  /rr  V  -  V  Y 

*1  (5a) 


IF.  r  ^n,  V  Y  as  of—  /w  V sin  cxP  =  /&?  Vy  — /7?1  Vof  (5b) 
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•  • 


My  =  —  I  $ 

C*1  /  (50 

These  inertia  terms  arise  when  an  aircraft  accelerates,  along  a 

curved  flight  path  as  shown  in  Fig.  4,  where  the  forces  and  moments  are 
expressed  as  d'Alembert's  inertia -resistance  reactions.  It  should  be  noted 
that  the  2  component  of  m$  namely,  the  contribution  -  m$  a,  is  small  ® 

compared  to  the  contribution  of  the  Z  component  of  centrifugal  force  mV?  and 
is  hereafter  omitted  from  Fg  ,  since  it  is  assumed  to  be  a  higher  order  term. 

No  similar  conclusions  can  iie  made  about  the  relative  size  of  the  X  com¬ 
ponents  of  these  forces  so  that  both  contributions  must  be  retained  in  F. 

Summing  all  the  forces  and  moments  as  stated  in  Eq.  1,  the 
longitudinal  dynamic  equations  of  motion  for  the  Flying  Platform  become: 
H-Force  Equation 


V 


(5a) 


(6b) 


EFx  =  -/>*■  Vy  —  o 

Tii fust  Equation 

£FZ  =  -T-c<fD  +  W+  'tjiVy  =  o 

Moment  Equation 

E  Yc/  n*-  H  •  ZH  +  TV  D  •  2D- 1 Y  £  -  o 

In  addition  to  the  small  perturbation  assumptions,  it  is  assumed  that 
the  changes  in  the  external  forces  and  moments  acting  on  the  Platform, 
because  of  the  small  departure  from  steady-state  motion,  will  depend 
entirely  on  the  displacement  and  disturbance  velocities  along  and  about 
the  aircraft  axis,  and  also  that  the  frequency  of  the  motions  is  low  so  that 
the  external  forces  and  moments  do  not  depend  on  the  accelerations  involved 


(6c) 
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Expanding  all  the  aerodynamic  forces  and  moments  of  the  former 
set  of  equations  as  a  Taylor  series  In  the  appropriate  variables,  but  neg¬ 
lecting  the  higher  order  terms ; 


rT 


1 2**+  — 


a«r 


aS 


D‘  *••• 


ffv  *  If* ^ f?*4^  •• 


and  letting: 


oS  -  +  A</ 


c<^=  ^  *  A(/ 


(9  «  +  AO 


S  -  A.  0 


V  =  V0  r  A  V 


y  = 


V  =  A  V 


e  -  of 


r  -  e~*r 


where  the  delta  variables  denote  tne  perturbations  and  a  ,  9  and  V 

000 

are  equal  to  zero,  since  it  is  assumed  that  the  motion  takes  place  about 
some  steady-state  trimmed  flight  condition  where  the  initial  rates  are 
zero. 

Also,  during  the  disturbed  motion  of  the  aircraft,  the  rotor  angular 
velocity  £2  is  assumed  constant,  so  that  Ajx  =  AV/£2R  . 
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Now  solving  the  longitudinal  dynamic  equations  (Eq.  6)  in  terms 
of  the  perturbations,  the  convenient  substitution  of  Eq.  6b  into  Eq.  6a 
may  be  made: 

/7n  /  i  «  7~  +  of  D  —  V/ 

(61 

so  that  Eq  6a  becomes: 

X  F  =  -/mV  -  o£  (T  ~w)  —  H  -  D~  WO  =  o  (9) 

X 

Substituting  Eq.  7  and  Eq.  8  into  Eq.  9  and  neglecting  higher  order  terms  , 
the  H-Force  equation  becomes: 


ZFX  -  -™v-H0-Da-wec-'F(ra-w  +  <F0D)-of'i^Lyt 

(10) 

For  the  initial  steady-state  condition  (subscript  o)  before  any 
acceleration  occurred,  Eq.  4  may  be  written  as: 


-H-D-we.  =  c 


ro  *  w  =  c 


U  la) 


(lib) 


+  £>■*=  o 

0  £) 


(11c) 


Substituting  Eq.  11a  and  Eq.  lib  in  Eq.  10  reduces  the  H-Force  equation 
to: 


#v» 
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0  ®  o>o^  ^ 


-  ~^e  —  4^  /ui  -  UA6 

<3<*r  a^  ^  / 


(12) 


2  2 

Nun-dimensionalizing  this  equation  by  dividing  Eq.  12  by  prrR  (i2R)  and 
treating  the  drag  term  as: 


r 


D  =  fcf  i/2 

p  7T$Z  (CiP)2  p  7Tfi>  z  (rh#f 


(13a) 


where: 


then: 


(1 3  Jo) 


(14) 


Also  introducing  the;  aerudynanuc  time  parameter: 


r  - 


/TTy 

prrfizruR 


and  using  the  notation: 


(15) 


)=  cLL-2 

}~  dp/r) 


(16) 


where  time  is  counted  in  the  dimensionless  time  ratio  t/r  ,  the  H-Force 


equation  in  the  usual  aerodynamic  form  is: 


(17) 


C_ 


+fc  +  c„  -d]  A©  *  0 

The  thrust  force  equation  may  be  developed  directly  from  Eq.  6b  using 
the  same  methods: 

^pn  V  Y  -  T-e£  D  +  W  =  O 

(6b) 


and  using  Eq.  15  and  Eq.  16  with  Eq  8: 

*  * 

dr  =  de-doY 
dG  *  d  C&B) 
cL°C  *  d  (&<d) 

the  thrust  force  equation  becomes: 


"c^c/ 
A 


o/i 


c 


*1 


fyL, 


V4C 


Co  +A 


Aid- 


AG  ~0  (18) 


For  the  moment  equation  (Eq.  6c)  non-dimensionalization  is  accomplished 

2  2  2 
by  using  prrR  (S2R)  R  ,  and  defining  1^=  mkY  or; 


and 


— 


ifi,T 


(19; 


where  k  is  the  radius  of  gyration  and  h^  is  the  non-dimensional 
moment  of  inertia  parameter. 
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The  moment  equation  is: 


%  *  r~Tu‘  jf-  ~  r-»  •  ~  t/L  C  Skjiui  +■  fc,  +  <Tr  .  *z 
A  a  ^  //?/««  Wy  c*r  *  1? 


A*  a 


.£r 


Sir  -&L*V 
V  rf  1 


£„  'd  -  C  •  %-<z-/>vd2]&e  =  o 

JU6)  die)  *  r  J 


Hence,  the  three  equations  of  motion  in  familiar  aerodynamic  form  are 
Eqs.  17,  13  and  20. 


H-Force  Equation 


- 'd  ('?-  y'rd  -  <//  ^ 


^  «  4*  M  V 


+  Fc  +  c 

j  W  ' 


"<z<c6)  ‘^1  &&  -  o 


Thrust  Equation 


Moment  Equation 


>  *  ^  ’  f  -  %  %  y<  C1  fj  A/*  +  & 


C  +cr.*r 

l*ls  ^  A 


r 

— ■  ’  _£T 


Acft^'d-r«  -itt.d-b  £2]a9  - 
L  ofcs;  /?  - 
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Explanation  of  the  Terms  of  the  Equations  of  Motion 

The  three  equations  of  motion  are  simultaneous,  homogeneous, 
ordinary,  differential  equations  with  constant  coefficients.  These  constant 
coefficients  are  expressed  by  stability  derivatives,  velocity  ratios,  a 
non-dimensional  inertia  parameter,  angle  measurements  and  drag  coeffi- 

c  c  c 

cients.  In  this  Report,  the  stability  derivatives  H  ,  H  ,  T 


a 


C  ,  c 

m  and  m 
a  a 


M  <*'  f 

are  evaluated  from  data  obtained  by  wind  tunnel 


experiment.  The  method  of  evaluating  the  stability  derivatives  and  drag 
terms  is  explained  in  the  Appendices. 

The  individual  terms,  expressed  in  non-dimensional  or  coefficient 
form,  in  the  equations  of  motion,  in  effect,  represent  the  following: 

Term 


1 .  CH  •  A  M 


2.  CR  -A  « 

a 

3.  CT  *  A  ii 


4.  CT  •  A* 
a 


represents  the  H-Force  due  to  the  change  Ap  in  the  volo- 

C 

city  ratio  where  H  is  the  rate  of  change  of  H-Force 

^  C 

coefficient  with  velocity  ratio.  For  the  case  of  H 

P 

H-Force  obtained  from  wind  tunnel  data  included  the 
measurement  of  all  the  aerodynamic  forces  of  the  rotor 
and  duct  assembly  acting  normal  to  ‘he  thrust  axis, 
represents  the  H-Force  due  to  a  change  A  a  in  the  angle 

Q 

of  attack  where  H  is  the  rate  of  change  of  H-Force 

a 

coefficient  with  angle  of  attack. 

represents  the  thrust  force  due  to  a  change  Ap  in  velo- 

Q 

city  ratio  where  T  is  the  rate  of  change  of  thrust  force 

A* 

coefficient  with  velocity  ratio. 

0 

For  the  case  of  T  ,  thrust  force  obtained  from  the 
wind  tunnel  data,  included  the  measurement  of  all  the 
aerodynamic  forces  of  the  rotor  and  duct  assembly  acting 
along  the  thrust  axis  , 

represents  the  thrust  force  due  to  the  change  Aa  in 

0 

angle  of  attack  where  T  is  the  rate  of  change  of  thrust 

a 


13. 


5 '  "oCD ' 


6  .  a  n  C_  1  A fj, 
oo  D 


7-IJtCD'A° 


8.  c  •  e*ti 

m 

/i 


9 .  C  • Ao 
m 

a 


force  coefficient  with  angle  of  attack, 
represents  the  drag  force  in  the  X  direction  contri¬ 
buted  by  the  pilot  due  to  the  change  Ap  in  the  velocity 

ratio.  The  term  C  '  equals  7— and  \x  is  the  initial, 

D  SR  o 

steady-state  velocity  ratio. 

(See  Appendix  I  for  additional  explanation  of  terms  5,6, 
and  7.) 

represents  the  drag  force  in  the  Z  direction,  contributed  by 
the  pilot,  due  to  the  change  Ap  in  the  velocity  ratio.  The 
initial,  steady-state  angle  of  attack  is  defined  aso^. 
represents  the  drag  force  in  the  Z  direction  contributed 
by  the  pilot  due  to  the  change  A  a  in  the  angle  of  attack. 

4 

represents  the  moment  contribution  measured  about  the 
quarter  chord  point,  due  to  the  change  Ap  in  the  velocity 

Q 

ratio  where  m  is  the  stability  derivative  denoting  the 
A* 

rate  of  change  of  pitching  moment  coefficient  about 

the  quarter  chord  point  with  velocity  ratio. 

The  stability  derivative  ^*m  was  evaluated  in  wind 

b 

tunnel  experiments  by  measuring  the  moment  created  by 
the  duct  and  rotor  only.  This  moment  was  arbitrarily 
transferred  to  the  quarter  chord  point  for  convenience 
in  accomplishing  additional  transfers  when  vertical  c.g. 
location  was  varied  in  order  to  study  the  effect  on  the 
stability  derivatives  .  Changes  in  vertical  c.g.  loca¬ 
tion  change  the  distances  e^,  x^,  and  e^. 
represents  the  moment  contribution  measured  about  the 
quarter  chord  point  due  to  the  change  A<*  in  the  angle  of 

Q 

attack  where  is  the  rate  of  change  of  pitching  moment 
coefficient  about  the  quarter  chord  point  with  angle  of 
attack.  This  derivative  was  evaluated  as  explained  in 
term  8 . 


* 


s.  at 


10. 


R 


20, 


C  Ap  represents  the  moment  contribution  about  the 

A 

^  center  of  gravity  due  to  the  thrust  force  T  •  Ap 

A* 


11 


R 


-CT  Aa 
a 


acting  at  the  moment  arm  x^. 

represents  the  moment  contribution  about  the  center 

Q 

of  gravity  due  to  the  thrust  force  T  •  A  a  acting  at  the 

a 


moment  arm  x 


r 


12. 

R  H 


■A p  represents  the  moment  contribution  about  the  center  of 

Q 

gravity  due  to  the  H-Force  H  •  Ap  acting  at  the 
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hL C 

R  H 


a 


moment  arm  z  . 

H 

ua  represents  the  moment  contribution  about  the  center  of 

Q 

gravity  due  to  the  H-Force  H  -A a  acting  at  the  moment 

a 


aim  z 


"D 


H' 


14.  — — p  C  ’A p  represents  the  moment  contribution  about  the  center  of 
R  o  D  » 

gravity  due  to  the  drag  force  pqCd'^a<  acting  parallel  to 
the  relative  wind.  The  perpendicular  distance  from  the 
center  of  gravity  to  the  line  of  action  of  this  drag  force 
is  the  moment  arm  z^ 

represents  the  weight  force  component  in  the  X  direction 

due  to  a  change  A0  in  attitude  angle  of  the  machine.  For 

the  component  C,.,9  in  the  X  direction  C,.r0  drops  out 
W  Wo 

from  Eq.  12  due  to  the  initial  condition  (Eq.  11a)  so  that: 


15 


cwAB 


or: 


C  6  =  C  9  +  C.  A0 
W  W  o  W 


C  8  =  C  A6 

W  W 


However,  in  the  thrust  equation,  the  weight  term  has  no 
perturbation  component  and  drops  out  entirely  from  Eq. 
6b  due  to  Eq.  lib. 


wi'fiiimw rnr'~inr  » imnrlr~Tfi-ii 


aiKM 
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16.  C. 


d(A9) 


■d(A@)  represents  the  additional  H-Force  due  to  the  linear 
velocity  caused  by  an  angular  pitch  rate  d(A0)  about  the 


17.  2H 


center  of  gravity. 

Q 

The  derivative  H  represents  the  rate  of 

d  (A0) 

change  of  H-Force  coefficient  with  pitching  velocity 
d  (A-jB) . 

**d^0''  rePresent:s  ^e  moment  contribution  about  the  center 


of  gravity  due  to  the  H-Force  h  , /a  ^.acting  at  the  moment 

d(A0) 


arm  z. 


This  term  can  be  thought  of  as  a  contribution  to 


18.  C 


‘d  (A  0) 


19.  d(A^) 


the  damping  in  pitch  of  the  aircraft. 
d(A0)  represents  the  moment  contribution  due  to  the  pitch 
rate  d(A0).  This  term  is  the  familiar  pitch  damping  term 
and  in  these  equations  refers  to  the  pitch  damping  cre¬ 
ated  by  the  rotor  and  duct.  Determination  of  this  value 
is  difficult  but  an  attempt  has  been  made  to  evaluate  the 
rotor  and  duct  contribution  to  pitch  damping  in  Appendix  II. 
represents  the  component  along  the  X  axis  of  the  inertia 
term  due  to  a  linear  acceleration  of  the  aircraft. 


20.  ^  d (A a)  -  fj.  d (A 0) • 
o  o 


•taken  together,  these  terms  represent  the 


component  along  the  Z  axis  of  the  centrifugal  force 
inertia  term  caused  by  the  rate  of  change  of  flight 


path  angle  7. 


21.  a  C_,  A p  +  a  Ct  Ao 
o  T  *"> 


-taken  together,  these  terms  represent  the 


component  along  the  X  axis  of  the  centrifugal  force  inertia 
term  caused  by  the  rate  of  change  of  flight  path  angle  7. 
The  form  of  this  term  is  due  to  the  substitution  performed 
in  Eq.  9  where  £q.  6b: 


/Tfi  V  i  =  T+  -  W 


(6b) 
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wad  substituted  into  the  H-Force  equation  as  a  matter 
of  convenience  in  solving  the  set  of  equations.  If 
Eq.  6a  had  been  solved  as  it  stands,  this  inertia  term 
would  have  been  in  a  lorin  similar  to  term  20  in  the 
thrust  equation.  However,  the  centrifugal  force  compo¬ 
nent  in  the  X  direction,  mVia,  of  Eq.  6a,  after  the  mathe¬ 
matical  substitution  and  operations,  takes  the  form  of 
term  2 1 : 

°^o  ^  ^/L  ^  -  Yf  oC 

This  indicates  that  this  bodv  axis  inertia  term  is  now 
represented  by  components  of  forces  parallel  to  the  rela¬ 
tive  wind.  All  terms  in  the  original  equations  represent 
the  summation  of  forces  and  moments  along  and  about  body 
axes  and,  in  this  case,  because  of  the  aforementioned 
substitution,  the  body  axis  inertia  force  is  represented 
by  components  of  forces  parallel  to  the  relative  wind. 
However,  these  still  represent  only  the  body  axis  inertia 
term  and  the  equations  are  still  consistent  with  the  body 
axis  frame  of  reference.  Changes  in  variables  and  sub¬ 
stitutions  among  the  set  of  Eq.  6  in  no  way  change  the 
validity  or  meaning  of  the  original  body  axes  equations 
of  motion  describing  the  dynamic  system. 

2 

.  lyr(Ae)  represents  the  pitching  moment  inertia  term  due  to  the 

distribution  of  the  mass  of  the  vehicle  where  h^  is  the 

non-dimensional  moment  of  inertia  parameter  measured 

2 

about  the  Y  axis.  The  term  d  (A0)  represents  the  per¬ 
turbation  angular  acceleration  about  the  Y  axis  measured 
with  respect  to  the  time  ratio  (t/r)  as  indicated  by  the 
symbol  (d) . 
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Solution  of  the  Equations  of  Motion 

In  order  to  solve  the  equations  of  motion  the  form  is  altered 
slightly  by  letting:  e  £v 


r 


Cl 


(21a) 


(21b) 


< 


C 


(21c) 


(2  Id) 


C,  -  C ,  +  *r  <?  r 


//  c. 


(2  i  e) 


(2  If) 


=  C 


<Zfr&) 


—  f// 

*  <£(*  e) 


(2  lg) 


where  the  bars  over  the  stability  derivatives  indicate  the  overall  or  total 
values  of  the  forces  and  moments  of  each  specific  derivative  acting  along 
or  about  the  axes  . 

From  the  above  notations  it  is  possible  to  see  that  changes  in 

C  C 

vertical  c.g.  location  will  affect  the  values  of  m  and  m  because  of 

M  o' 


e 


changes  In  zTT,  z  ,  and  xm.  The  term  m  , ,  A  will  not  change  radically 
c  H  d  T  d(A0) 

since  is  a  rather  small  quantity. 

Also,  the  moment  of  inertia  about  the  Y  axis  L^,  is  held  constant 
for  this  dynamic  analysis. 

Using  this  notation,  Eqs.  17,  18  and  20  take  the  form: 

H-Force  Equation: 


-'jV' 


f  c#  •  d.\&e 

d(*e)  J 


(22a) 


T-Force  Equation 


[%]*/*  ~  fa cL]*Q  *  0 


(22b) 


Moment  Equation 


P>]**  2  -  ° 

The  equations  of  motion  are  three  simultaneous  differential  equa¬ 
tions  with  constant  coefficients,  the  solutions  of  which  may  be  obtained 
by  assuming  a  solution  in  the  form: 

A%  ,  aW 

aG^aQ^ 


(22  c) 


A  'T 


and,* 


d(&6 )-  AfyJe 


*  % 


2  .2^-^r 

d  (&$)=  aGjA  e 


where  X  is  a  real  or  complex  constant  of  equal  value  for  each  variation, 
and  where  A^,  ,  and  A0  ,  are  also  real  or  complex  constants. 
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Substituting  Eq.  23  and  Lq.  24  into  £q.  22,  and  dividing  out 
X  t/j 

the  common  term  e  '  ,  reduces  Eq  22  to  three  algebraic  equations 

in  the  unknown  X  and  the  terms  A/u^ ,  Aa^  and  A0^ ; 


[%'J] y1'  +[c**<w-ik9'" 


G  (25a) 


_%]aa  -  \ao^]a&'  =  0 


(25  b) 


r  1  c 

H,  aa  +  ■ 


^  +  C, 


»  1 

d(t ie) 


hy  A  A9=  O 


(25  c) 


It  is  a  property  of  this  type  of  equation  that  there  can  be  non- zero 
values  of  the  unknowns  only  when  the  determinant  of  the  coefficients  is  ° 


equal  to  zero: 


G +  £  '  ^ 

"  <?A6) 


CrAA°*  A  ^  =  o 


C, 

/9yu. 


C 


C?  “  Ay  2 
<*£e)  Y 


Expanding  this  stability  determinant;  the  characteristic  equation 


in  X  is  of  the  form: 


A\*+BX*+cf  +&X+E  —  o 


vam  mm. 


a-  t i 
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where : 


/?  ••=  Ay/** 


(28a) 


B- 
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y° 

C  - 

-  <?r 

o*r 

?  ~m.b  C  <L 

<&»)  ^  ^  <?£#) 

£ 

9y 

c„ 

<f(*8) 

£ 3  ~/t> 

c  ^ 

y  p* 

> 

'"’ike) 

+/*• 

-  c  c  c 

+/to 

£ 

c 

-A> 

CL 

y 

+  crj  c„ 

5 

6) 

£= 

(Gy 

-  C  r  ' 

1 

(28b) 


(28c) 


(28d) 


(28e) 


It  is  convenient  at  this  point  to  consider  the  hovering  Flying 
Platform  as  a  special  case  of  the  former  development. 

The  common  assumption,  in  hover  analysis,  is  the  exclusion  of 
the  vertical  degree  of  freedom  since,  in  hovering,  the  horizontal  motion  and 

c  c 

tilting  do  not  affect  the  thrust.  The  derivatives  T  and  T  are 

H  a 

are  zero.  Hence,  the  thrust  will  be  essentially  constant  and  there  will 


■r  MMwMiWg  mi  aiwiUjL.  ijfiir  y  wwiii  KwiMXh 


27 


be  no  appreciable  vertical  motion  of  the  machine,  principal  interest 
is  centered  cn  changes  of  velocity  ratio  and  attitude  angle  and  not  on 
altitude  change.  This  can  be  accomplished  mathematically  by  altering 
the  original  equations  of  motion  so  that  the  thrust  force  equation  is  el i,n- 
mated  entirely.  The  derivatives  CT/i  and  CTff  are  zero  and  the  alpha 
te-ms  in  the  H-Force  and  moment  equations  are  dropped. 

The  algebraic  equations  are: 


and  thedeterminant  of  coefficients  is: 


<?(*&) 

/  ( Gp?  *  ^  -  /)^  J 

<?(&&)  y 


(29) 


A9,  =  (9 


c 


=  O 


(30) 


Expanding  this  determinant  yields  the  third  order  characteristic 
equation: 

*  £J*  +  CJ  +  £>  =  o 


(31) 


where: 


(32a) 
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(32b) 


(32  c) 


(32d) 


The  characteristic  equations  obtained  for  forward  flight  and 
hovering  are  often  referred  to  as  the  stability  quartic  and  cubic.  The 
roots  of  these  equations  are  the  values  of  X  that  determine  the  final  solu¬ 
tion. 

Then,  for  example,  the  response  of  an  aircraft  in  pitch  AG  to  a 
disturbance  may  be  found  by  a  solution  of  the  following  form: 

Ifr  XU  X/Sr 

A6  =  a6,S  i-A62c  ,A&se  +ABye 


Determination  of  time  histories  of  the  variables  by  actual 
solutions  of  this  form  is  quite  lengthy  and  tedious  if  numerous  deriv¬ 
ative  variations  are  made . 

However,  analog  computers  can  be  used  to  solve  the  equations  of 
motion  in  order  to  obtain  the  time  histories  associated  with  arbitrary 
changes  in  the  stability  derivatives.  Typical  time  histories  are  shown 
in  Figs .  5  and  6 . 
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Investigation  of  the  Roots  of  the  Characteristic  Equation 

Frequently,  the  important  aspect  in  the  study  of  a  dynamic  system 
is  the  determination  of  the  characteristic  modes  of  motion  of  the  aircraft 
after  a  disturbance  from  steady  flight,  by  an  investigation  of  the  four 
roots  of  the  stability  quartio. 

Methods  of  predicting  the  character  of  the  motion  from  the  roots 
or  the  coefficients  are  covered  in  various  texts.  (Ref.  5). 

It  is  of  interest  to  know  whether  the  motion  is  oscillatory  or 
aperiodic.  If  oscillatory,  it  is  necessary  to  determine  the  period  and 
datnping.  When  aperiodic,  it  is  important  to  determine  the  rate  of 
divergence  or  convergence. 

This  information  may  be  obtained  by  examining  the  roots  of  the 
stability  quartic. 

When  the  roots  are  positive  real  numbers,  the  motion  is  an  aperi¬ 
odic  divergence;  negative  real  numbers,  aperiodic  convergence;  and 
complex  pairs,  an  oscillatory  motion.  In  the  case  of  a  complex  pair  of 
roots,  the  oscillation  is  undamped  if  the  real  part  is  positive,  and  damped 
if  the  real  part  is  negative. 

The  roots  of  the  quartic  may  be  obtained  by  any  of  the  numerous 
methods  outlined  in  many  texts . 

In  general,  if  the  roots  are  both  real  and  complex  as; 


the  period  of  the  oscillatory  pair  may  be  found  by; 


p  —  ^  r* 

l V 


[seconds] 


(35) 
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and  the  damping  may  be  found  as  the  time  (T)  required  for  the  transient 
response  to  grow  to  double  amplitude  or  to  decay  or  damp  to  half 
amplitude: 

7^=  —  — 7"  =■  93  pf  r 5 £cc/?c/z] 

Cb  Ob  L  (36) 

Taking  the  center  of  gravity  at  the  original  location,  (.42c  above 
the  top  of  the  duct)  the  roots  of  the  quartic  for  the  different  ducts  at 
different  speeds  are  computed .  The  damping  and  period  are  shown  for 
each  of  the  roots.  See  Tables  I  and  II. 

Investigation  of  the  Coefficients  of  the  Characteristic  Equation 

Considerable  information  concerning  the  character  of  the  motion 
of  the  aircraft  can  be  obtained  from  the  coefficients  of  the  stability  quartic 
without  solving  for  the  roots . 

If  all  the  coefficients  of: 

=  O  (37) 

are  positive,  there  can  be  no  positive  real  root  and  there  is  no  possibility 
of  a  pure  divergence.  If  one  of  the  coefficients  changes  sign  (becomes 
negative)  there  cap  be  either  an  increasing  oscillation  or  a  pure  diver¬ 
gence  in  one  of  the  modes  . 

The  C  coefficient  is  of  special  importance  since  it  represents 
approximately  the  maneuver  margin  of  the  aircraft  (Ref.  6).  If  the  C 
coefficient  changes  sign  (become^  negative)  it  is  an  indication  of  a 
serious  instability.  There  will  probably  be  a  rapid  divergence  or  unstable 
oscillation. 
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The  £  coefficient,  where: 

r  3C.  5Cr  3C,  d£r 

^  L  3A  W  3c*r  (38) 

is  of  importance  for  static  stability.  The  static  stability  is  positive  if  E 
is  greater  than  zero.  If  E  is  equal  to  zero,  one  of  the  roots  is  zero  and 
one  of  the  modes  can  continue  unchanged  indefinitely.  Usually  E  is  the 
first  coefficient  to  change  sign,  and  therefore  the  first  indicator  of  a 
divergence  or  positive  real  root.  In  some  instances,  it  is  possible  for  C 
kto  change  sign  before  E,  in  which  case  there  will  be  a  dynamic  instability 

r 

occurring  with  a  positive  static  stability. 

Regardless  of  the  signs  of  the  coefficients,  it  is  possible  to  deter¬ 
mine  the  number  of  roots  having  positive  real  parts  by  using  Routh's 
criterion.  This  method  requires  only  simple  calculations  utilizing  the 
coefficients  of  the  stability  quartic  in  setting  up  the  familiar  Routh's 
array.  The  number  of  roots  having  positive  real  parts  is  equal  to  the 
number  of  sign  changes  occurring  in  the  first  column  of  the  array. 

Also,  from  Routh's  array,  it  is  possible  to  obtain  the  combination 
of  coefficients  known  as  Routh's  discriminant: 

R  =  BCD-/)DZ-BS£  (39) 

where  A,  B,  G,  D  and  E  are  the  coefficients  of  the  quartic.  If  all  the 
coefficients  are  positive,  it  is  possible  to  make  the  followino  statements 
concerning  the  stability  of  the  system. 

If  Routh's  discriminant  is  positive,  there  is  no  possibility  of  the 
real  part  of  any  complex  pair  being  positive.  There  will  be  no  undamped 
oscillations . 
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If  Routh'e  discriminant  is  negative,  there  is  a  complex  pair  of 
roots  with  a  positive  real  part.  There  is  an  undamped  oscillation. 

If  Routh's  discriminant  is  equal  to  zero,  there  is  a  neutrally  “damped 
oscillation  implying  a  complex  pair  of  roots  with  a  zero  real  part. 

The  criterion,  (R  =  0),  then  indicates  &  boundary  between  stable 
and  unstable  oscillations. 

This  condition,  (R  =  0),  developed  from  the  coefficients  of  the  sta¬ 
bility  quartic,  may  be  investigated  further  by, assuming  a  solution  of  the 
form  for  neutral  oscillations  (X  =  iw).  Substituting  this  solution  with  the 
zero  real  part  into  the  characteristic  equation  gives: 

A)*+B13  +  CA%£X  +£  =  O 

/) / CO*  t £i3603  t-  Ci*U)2  +  /)/&  /-  S  =  O  (40) 
&cu3i  -  Coj*  +Dtui  t-  £  —  o 


Forming  two  equations  of  the  real  and  imaginary  parts : 

/?U>  *  —  C  U)  *  7*"  -  o 

-  £oj*  +  Doi  -  o 
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orj 


Substituting  Eq.  43  in  Eq.  41  results  in: 


(45) 


SC/P  -AD2-82£  =  0-8 


(46) 


which  is  the  form  of  Routh's  discriminant. 

With  this  analysis,  one  can  obtain  an  appreciation  of  the  require¬ 
ments  imposed  on  the  coefficients  A,  B,  C,  D,  E,  so  that  the  solution 
R  =  0  is  a  stability  boundary  between  stable  and  unstable  oscillations. 

Since  a  solution  of  the  form  for  neutral  oscillations  is  assumed, 

(X  =  iu),  inspection  of  Eq.-43  reveals  that  *  will  be  a  real  number,  e.g.  , 
u  =  a,  if  D  and  B  are  both  positive,  and  the  solution  X  =  ia  will 
indicate  a  neutral  oscillation.  However,  if  D  or  B  changes  sign  so 
that  u  =  ia  in  Eq.  43,  then  the  solution  will  be  X  =  -  a,  or  a  real  root. 
Therefore,  when  the  D  coefficient  changes  sign,  (becomes  negative), 
the  solution  of  Routh's  discriminant  equal  to  zero  no  longer  defines  an 
oscillation  boundary. 

Boundaries  of  the  Stability  Quartic 

In_the  analog  computer  program,  it  was  found  that  the  derivatives 

c  c 

m  and  m  had  a  principal  effect  on  the  mode  of  motion  of  the  Flying 

hi  a 

Platform. 

The. Figs.  10A,  11A,  13Aand  14A  entitled  "Boundaries  of  the 

C  C 

Stability  Quartic"  are  drawn  with  m  and  m  as  the  ordinate  and 

ji  a 

abscissa. 
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Included  bn  this  diagram  are  various  boundaries  of  the  stability 

.  Q  Q 

quartic  plotted  as  functions  of  m  and  m  .  Values  of  the  total  stability 
derivatives  were  used  as  calculated  from  the  wind  tunnel  data  for  each 
speed  and  duct  configuration.  The  theoretical  value  of  pitch  damping 
was  used  and  kept  constant  for  these  diagrams. 

The  coefficients  o i  the  stability  quartic,  G,  D  and  E  are  plotted 

c  c 

cis  functions  of  m  and  m  and  appear  as  the  straight  lines  labeled 

a 

C  =  0,  D  =  0  and  E  =  0.  The  lines  so  labeled  represent  the  boundaries  ajt 
which  each  coefficient  changes  sign  from  positive  to  negative  for  the 

c  c 

appropriate  values  of  m  and  m  . 

li  a 

The  coefficients  are  negative  on  the  hatched  side  of  the  line.  For 
example,  near  the  origin  in  the  right  half  plane  the  vertical  line  labeled 
C  =  0  is  negative  on  the  right  side  of  the  line. 

The  meaning  of  the  traversing  of  these  lines  (coefficients  changing 
sign  from  positive  to  negative)  has  been  discussed  previously  and,  for 
this  section,  can  be  summarized  as: 

C<  0  -  Serious  instability  associated  with  negative 
maneuver  margin. 

D  <  0  -  The  solution  R  -  0  no  longer  has  meaning  as 
an  oscillation  boundary. 

E"<  0  -Instability  associated  with  a  negative  static 

margin.  Aircraft  statically  unstable.  Dynamic 
instability  present  in  the  form  of  a  divergence 

due  to  static  instability. 

C  C 

The  coefficients  A  and  B,  not  functions  of  m  and  m  and 

M  « 

generally  positive,  do  not  appear  on  these  diagrams. 

The  line  R  =  0,  representing  Routh's  discriminant,  is  plotted 
to  denote  the  boundary  between  stable  and  unstable  oscillations. 
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c  c 

This  line  represents  the  loci  of  the  values  of  m  and  m  for 

a 

neutral  oscillator/  motion. 

On  die  left  Side  of  the  line  R  is  positive  and  the  oscillations  will 
be  stable  .  On  the  right  or  hatched  side  R  is  negative  and  the  oscillations 
will  be  unstable.  The  line  R  =  0  does  not  exist  past  the  line  D  =  0,  since 
its  meaning  as  an  oscillation  boundary  is  ended  whenever  the  D  coefficient 
becomes  negative. 

In  one  case,  (Fig.  13  A),  the  line  R  =  0  has  essentially  three 
branches,  all  representing  a  neutral  oscillation  boundary,  and  all  located 
in  the  region  where  the  D  coefficient  is  positive.  In  the  top  half  plane, 
the  two  branches  of  R  =  0  form  a  parabola  ineide  of  which  there  is  an 
unstable  oscillation  (at  least  two  of  the  roots  must  form  a  complex 
pair  with  a  positive  real  part).  In  the  bottom  half  plane,  only  one 
branch  of  the  parabola  appears  since  the  other  is  nonexistent  when  it 
crosses  the  D  =  0  line  into  the  region  where  the  D  coefficient  is  negative. 
Note  that  the  line  E  =  0  passes  through  the  middle  of  the  R  =  0  parabola  in 
the  top  half  plane . 

c  c 

Whenever  the  values  of  m  and  m  locate  a  point  on  the  negative 

a 

side  of  the  E  =  0  line,  at  least  one  of  the  roots  of  the  quartic  must  be 
positive  end  real,  and  an  instability  exists  which  takes  the  form  of  a 
divergence . 

Then,  In  the  region  between  the  E  =  0  line  and  the  right  hand  branch 
of  the  R  =  0  line  in  the  first  quadrant  (Fig.  13A),  at  least  three  of  the  roots 
of  ttfe  quartic  are  accounted  for;  a  complex  pair  with  a  positive  real  part 
and  a  positive  real  root.  The  other  root  in  this  region  is  a  convergence 
as  seen  in  the  root  It  cus  plot  Fig.  13-1. 

The  other  two  branches  of  the  R  =  0  line  located  in  the  second  and 
third  quadrants  represent  oscillation  boundaries,  so  that  in  the  region 
between  them,  R  is  greater  than  zero  and  the  entire  area  represents  a 
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completely  stable  region  of  only  convergences  and  stable  oscillations. 

In  Fig.  14A,  the  R  =  0  line  is  terminated  when  it  reaches  the  E  =  0 

Q 

line  since  it  passes  into  a  region  where  m^  and  E  are  always  negative, 
so  that  the  predominating  mode  of  motion  is  a  divergence  that  strongly 
masks  all  other  modes  of  motion. 

The  line  labeled  RLBP  represents  a  locus  of  points  obtained  from 
root  locus  diagrams .  The  letters  RLBP  are  an  abbreviation  for  "Root  Locus 
Break  Points".  The  break  point  in  this  case  represents  the  point  in  the 
right  half  plane  of  the  root  locus  diagram  where  a  complex  pair  of  roots 
joins  the  positive  real  axis,  thereby  terminating  the  unstable  oscillation 
and  becoming  two  aperiodic  divergences. 

In  other  words,  moving  to  the  right  from  the  stable  region,  across  the 
R  =  0  line  ,  and  passing  into  a  region  of  unstable  oscillations,  at  least 
two  of  the  roots  of  the  quartic  are  a  complex  pair  with  positive  real 
parts  (Figs.  10A,  11  A,  and  14A).  Continuing  to  move  to  the  right. into 
the  first  quadrant  and  crossing  the  RLBP  line,  this  Complex  pair  becomes 
two  real  positive  roots,  ending  the  unstable  oscillation  and  now  'epre- 

■j. 

sentinq  two  divergent  roots .  To  the  immediate  right  or  hatched  side  of 
the  RLBP  line,  at  14a  st  two  of  the  roots  of  the  quartic  are  positive  real 
roots . 

The  line  RLBP  does  not  appear  in  Fig.  13A,  because  the  parabola 
R  =  0  already  defines  a  large  region  where  unstable  oscillations  start 
and  stop. 

Q  Q 

However,  for  extremely  large  values  of  m  and  m  beyond  the 

a 

values  of  this  diagram  (Fig.  13A),  there  is  a  region  where  the  complex 
roots  become  positive  real  roots  for  a  time.  (See  root  locus  Fig.  13-I-e, 

c  c 

points  180  and  181.)  The  value  of  m  is  equal  to  0.006  and  m  is 

a  [i 

approximately  equal  to  0.1  at  these  points,  which  &re  beyond  the  pertinent 
C  C 

range  of  m  and  m  for  this  particular  case. 

a 
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The  various  center  of  gravity  locations  have  been  plotted  on  these 

c  c 

diagrams  according  to  their  respective  values  of  m  and  m  as 

/ L  O' 

listed  in  Table  III  .  The  broken  line  between  these  points  represents  the 
approximate  line  of  travel  for  the  c.g.  at  intermediate  values. 

Explanation  of  Stability  Diagrams 

The  Figs.  10B,  11B,  13B  and  i4B  entitled  “Stability  Diagrams " 
are  a  greatly  Simplified  form  of  the  previous  boundary  diagrams.  The 
only  portions  of  the  boundaries  that  have  been  retained  are  those  which 
divide  the  diagram  into  the  pertinent  regions  of  stability. 

The  unshaded  area  represents  a  region  of  complete  dynamic  sta¬ 
bility  .  Only  roots  representing  convergences  and  stable  oscillations 
exist  in  this  region. 

The  shaded  areas  represent  regions  of  dynamic  instability.  Roots 
existing  in  this  region  represent  divergences  and/or  unstable  oscillations. 

There  are  always  four  roots  to  contend  with  in  every  region.  In 
the  stable  region  all  four  roots  must  be  stable  roots.  In  the  unstable 
region,  a  minimum  of  one  unstable  root  must  always  exist  while  the  others 
may  be  stable. 

Actually,  the  unstable  roots  are  the  predominant  modes  of  motion  In 
the  unstable  regions  and  mask  the  stable  modes  present. 

One  shaded  unstable  region  is  labeled  "unstable  oscillations" 
and  indicates  that  at  least  two  of  the  four  roots  in  this  region  must  be 
a  complex  pair  with  a  positive  real  part.  The  other  shaded  area  indicates 
the  unstable  region  where  a  minimum  of  one  unstable  root  must  exist. 

Some  of  the  interesting  modes  of  motions  which  prevail  in  the 
various  regions  are  labeled  in  these  diagrams. 


m 


1  o 
Oc 

Again,  the  various  center  of  gravity  locations  have  been  plotted 

c  c 

according  to  their  respective  values  of  m  and  m  as  listed  in  Table 

a 

III.  The  broken  line  represents  the  approximate  line  of  travel  of  the  center 
of  gravity  at  intermediate  values . 

C 

For  the  case  of  hovering,  u.  and  m  equal  zero.  The  diagram 

C  °  a 

is  plotted  only  for  m  and  appears  in  the  form  shown  in  Fig.  8. 

The  study  of  these  diagrams,  together  with  the  associated  root  locus 

diagrams,  will  present  quite  clearly  and  concisely  the  precise  modes  of 

motion  to  be  expected  at  different  speeds  for  different  center  of  gravity 

locations  with  the  different  duct  configurations. 

# 

Explanation  of  the  Root  Locus  Diagrams 

The  stability  diagrams  were  prepared  in  order  to  show  concisely 

c  c 

and  clearly  the  role  of  m  ,  m  ,  vertical  center  of  gravity  travel,  speed 

/i.  a 

changes,  and  different  duct  designs,  in  affecting  the  modes  of  motion 
of  the  Flying  Platform. 

The  diagrams  exhibit  rathter  precisely  the  mode  of  motion  associated 
with  locations  in  the  various  regions . 

Although  the  determination  of  the  stability  boundaries  is  quite 
lengthy,  they  may  be  calculated  by  various  methods. 

In  the  formulation  of  the  stability  diagrams  for  this  Report,  use 
was  made  of  the  root  locus  method  of  displaying  dynamic  responses. 

With  this  method,  it  is  possible  to  exhibit  the  properties  of  a  linear 
system  by  plotting  the  loci  of  the  roots  of  the  characteristic  equation  in 
which  the  values  of  selected  derivatives  have  been  varied. 

A  complete  treatment  on  the  utilization  of  this  method  is  avail¬ 
able  in  texts  by  Evans  (Ref.  7)  or  Savant  (Ref.  8). 
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Since  primary  interest  is  centered  on  varying  the  derivatives  m 

C  1 

and  m  ,  the  characteristic  equation  may  be  written  in  the  form: 
a 


+  C  £ 

Hz 


C  C 

3  ^ 


)J 


where  the  quantities  enclosed  in  the  parenthesis  are  equal  to  the  coeffi¬ 
cients  C,  D  and  E.  From  this  expression  it  is  possible  to  calculate  the 
polas»nd  the  zeros.  Since  vertical  center  of  gravity  travel  seems  to  have 

c 

the  primary  effect  of  changing  m  ,  the  root  locus  diagrams  were  plotted 
C  ^ 

so  that  m  varied  from  zero  to  positive  or  negative  infinity.  The  poles 

M  C  C 

were  calculated  for  specific  values  of  m  with  m  equal  to  zero. 

a  il 

Solving  the  characteristic  equation  in  this  manner  gives  the  four  roots  of 
fhequartic  p  ,  p  ,  p  ,  p  .  The  zeros  of  the  equation  are  the  increments 

Q  X  Z  0  T 

of  m  ,  namely: 


C3 


4 


=  c 


IV 


(48) 


Q 

Solving  the  above,  the  zeros  for  m  are  obtained  and  denoted  by  z. 

C  ^ 

Similarly,  the  zeros  for  m  denoted  by  Z,  may  be  obtained  from: 

a 


4  - 

d,  =  -  Cr  cM 


(49) 


l 
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i 
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Using  the  former  expressions,  the  root  locus  equation  may  be  writ- 

c  c 

ten  for  a  specific  value  oi  m  with  m  varying  from  zero  to  positive 

a  p. 

or  negative  infinity  as: 


where  the  p's  represent  the  poles  and  the  z's  the  zeros  of  m  . 

C  P  Q 

For  the  case  where  the  poles  (?)  are  calculated  for  m  =  v,/m 

C  M  a 

and  the  zeros  (Z)  for  m  are  utilized,  the  expression  takes  the  form: 


=  0, 


4  m  -  p3 )  -  c  -z,)u-*yo  <5  u 


c  c 

in  which  m  varies  from  positive  to  negative  Infinity  while  m  equals 
zero. 

A  group  of  root  locus  diagrams  are  plotted  so  that  essentially  all 

c  c 

values  of  m  and  m  are  covered  for  a  particular  duct  at  a  selected 

pa 

speed. 

Associated  with  each  stability  diagram  there  are  five  root  locus 
sets.  The  first  set  investigates  the  modes  of  motion  encountered  when 
moving  along  the  ordinate  or  abscissa  of  the  stability  diagram.  The 
other  four  sets,  labeled  I,  II,  III  and  IV,  investigate  the  modes  of  motion 
present  in  the  four  quadrants  of  the  stability  diagram.  The  value  of 

Q 

m^  (Eq.  50)  in  these  diagrams  varies  from  zero  to  positive  or  negative 
infinity. 

The  root  locus  sets  are  plotted  for  both  ducts  at  pi  equal  to  0.05 
and  0.10. 
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Included  in  each  figure  used  to  investigate  the  modes  of  motion  in 
each  quadrant  (I,_II,  III,  IV)  ,  one  root  iocu3  is  plotted  with  the  zeros  of 

c  c 

m  instead  of  m  .  This  plot  is  divided  from  the  rest  by  a  double  lin^ 
a  At 

border  as  in  Fig.  10-I-f.  In  each  of  these  particular  double  border 

c  c 

diagrams,  m  is  equal  to  zero  and  fn  varies  from  zero  to  infinity.  This 
[L  a 

then  shows  the  movement  of  the  initial  positions  of  the  poles  for  each  queri¬ 
ed 

rant  as  the  value  of  m  is  increased  from  zero  to  infinity.  In  this  way  it 

«  C 

is  possible  to  sketch  the  root  loci  for  all  higher  values  of  m  beyond 

a  C 

those  calculated  for  each  quadrant.  In  this  respect,  the  values  of  m 

a 

in  each  quadrant  are  usually  large  enough  so  that  the  basic  appearance 
of  the  root  loci  will  not  change  for  values  beyond  the  highest  shown  in 
each  set. 

All  root  loci  drawn  in  this  report  were  accurately  calculated  and 
plotted  to  the  scale  shown  in  each  figure. 

The  root  locus  sei  for  =  0  is  also  plotted  in  Fig.  7.  For  both 

Q 

ducts  in  the  hovering  case  m  equals  zero;  The  stability  cubics  are 

a 

written  as; 


(31) 


The  root  locus  equation,  after  dropping  a  negligible  "zero"  due  to 

Q 

the  small  value  of  H,  ,  *  ,  is; 

d  ( AQ) 


(53) 


MR 
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where  C^.  is  known  and  m^  is  varied  from  zero  to  positive  or  negative 

Q 

infinity.  The  poies  are  calculated  for  m  equal  to  zero.  The  diagrams 
for  the  hovering  case  have  three  poles  and  no  zeros . 

The  root  locus  method  may  also  be  utilized  conveniently  for  the 
investigation  of  the  effects  of  artificial  stabilization  on  the  transient  re-- 
spOnseof  the  Flying  Platform.  Artificial  stabilization  may  be  provided 
by  a  rate  sensitive  device  such  as  the  gyro-bar  with  aerodynamic  or 
mechanical  damping. 

The  addition  of  thifl  type  of  device  to  the  dynamic  system  may  be 
accomplished  mathematically  by  ir  eluding  a  rate  sensitive  transfer  function 
with  a  phase  lag  of  the  type  --  ---  •  —  to  the  pitch  damping  term  in  the 

moment  equation. 

For  hovering,  the  determinant  of  coefficients  in  algebraic  form 

* 

becomes : 


Expanding  this  determinant: 


0  =  Ar/l  +  (by  \)  ^CyjC  +A1_  (x+C  \  , 

^  s'  l+Tr  (55) 


and  arranging  the  poles  and  zeros  as: 


where  T  is  the  time  constant  for  the  phase  lag  ,  the  root  locus  equation  for 
the  hovering  case  may  be  written  as: 


4-  ( Yh )  =  °  <5  7> 


where  the  gain  of  the  feedback  loop,  K,  of  the  stabilizing  device  varies 
from  zero  to  infinity. 

The  determinant  of  coefficients  for  the  addition  of  artificial  stabili¬ 
zation  in  forward  flight  is: 


c 

*eC 

c  +  c 

■A 

d  Y*°J' 

-/*•* 

=  0 

(58) 

c 

d 

AT. A 
UTX 

and  expanding  the  determinant: 


6=/tl  ¥+  31  \  C/+  HI  +£  + 


UT1 


(59) 


where  the  first  five  terms  represent  the  original  characteristic  equation 
and  T  is  again  the  time  constant  for  the  phase  lag.  Separating  the  poles 
and  zeros  in  root  locus  form  gives: 


/*o  Ay  *  £ (*)(!-  ?t)U-  =  ° 


(60) 


»»•« 
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where  K,  the  gain,  varies  from  zero  to  infinity. 

The  root  locus  diagrams  for  the  artificially  stabilized  machine  in 
hovering  and  high  forward  speed  are  plotfed  in  Figs,  9,  12  and  15. 

These  diagrams ,  together  with  all  other  pertinent  root  locus 
diagrams,  are  fully  explained  in  the  DISCUSSION. 


£> 


Effect  of  Vertical  Center  of  Gravity  Movement  on  the  Overall  Stability 
Derivatives 

The  overall  stability  derivatives  are  expressed  by  £q.  21  as: 

(2  If) 

(21e) 


C, 


C'  +  Zr  c 


4  4 


- 


A 


4 


The  distances  x^,,  and  zQ  represent  the  moment  arms  from 
the  total  center  of  gravity  to  the  respective  forces,  ^m  and 

#1  a  q 

lepresent  the  stability  derivatives  about  the  quarter  chord,  whereas,  m 

C  ^ 

and  m  are  functions  of  the  distances  x.,,  z„  and  z^ ,  When  the  total 
q-  T  K  D 

center  of  gravity  of  the  aircraft  is  moved  vertically  these  distances  will 
vary  so  that  considerable  change*  in  the  overall  value  of  the  derivatives, 

c 

especially  m  ,  can  be  obtained. 

M 

In  the  original  c.g.  location,  the  total  center  of  gravity  is  0.42 
of  the  chord  length  (0.42c)  above  the  top  edge  of  the  duct.  For  the  case 
of  the  long  chord  duct,  the  total  center  of  gravity  of  the  aircraft  is  located 
approximately  nine  inches  above  the  top  edge  of  the  duct.  Similarly,  the 
center  of  gravity  of  the  short  chord  duct  is  also  assumed  to  be  at  the  0.42c 
location  so  that  it  is  approximately  five  inches  above  the  top  of  the  duct. 
These  positions  are  designated  as  the  original  center  of  gravity  positions 
from  which  the  vertical  variations  of  c.g.  are  measured.  The  pilot's 


sta  id-on  platform  always  remains  at  this  original  position,  i.e.,  0.42c 
above  the  top  of  the  long  and  short  chord  duct. 

The  actual  scale  profiles  of  the  long  and  short  chord  duct  models 
are  shown  in  Fig.  2 . 

c  c 

The  values  of  the  stability  derivatives  m  and  m  for  the  different 

ii  a 

center  of  gravity  locations  shown  in  Fig.  ]6  are  presented  in  Table  III. 


c.g.  0. 5  ft.  original 

down  from  location 

original  of  c.g. 


c.g.  0.5  ft.  up 
from  original 
location 


c.g.  1 . 0  ft.  up 
from  original 
location 


Fig.  16  CENTER  OF  C  7ITY  LOCATIONS 
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DISCUSSION 


General 


1.  Similarity  of  Stability  Derivatives 

From  the  information  contained  in  Table  III,  concerning  the  over¬ 
all  value  of  the  stability  derivatives,  it  is  apparent  that  there  is  some  sim¬ 
ilarity  or  overlapping  of  values  for  the  two  ducts  for  different  speeds  at 

Q 

different  c.g,  locations  .  For  example,  the  values  of  m^  for  the  long 
chord  duct  at  p  equal  to  0.05  with  the  c.g.  raised  0.5  feet  is  equal  to 

Q 

0.0136  while  m^  for  the  short  chord  duct  at  n  equal  to  0.05  with  the  c.g. 

at  the  original  location  is  approximately  the  same  or  equal  to  0.0121. 

A  small  adjustment  in  c.g.  location  of  one  duct  configuration  would  make  the 

derivative  equal  that  of  the  other  configuration.  Similarly,  the  value  of 
C 

m^  equal  co  0.0032  also  matches  up  for  both  ducts  when  /i  equals  0.10 
at  the  highest  c.g.  location. 

Therefore,  at  a  specific  speed  in  a  certain  speed  range,  similar 

m 

values  of  either  stability  derivative  may  be  obtained  for  both  duct  con¬ 
figurations  by  adjusting  the  vertical  center  of  gravity  location. 

2 .  Attitude  Angle  for  Steady  Forward  Flight 

As  seen  in  Fig.  18,  the  attitude  angle  for  steady-state  equilibrium 
in  forward  flight  is  plotted  versus  /*.  For  convenience,  Fig.  17  is  plotted 
directly  above  Fig.  18  so  that  conversions  may  oe  easily  accomplished 
between  and  V  in  feet  per  second  or  miles  per  hour. 

From  Fig.  18,  it  is  evident  that  the  long  chord  duct  requires  larger 
tilt  angles  than  the  short  chord  duct.  This  can  be  explained  by  the  fact  that 
the  parasite  and  momentum  drag  of  the  long  chord  vehicle  is  probably 
greater  than  that  of  the  short  chord  duct. 

For  a  long  chord  duct,  the  flow  tends  to  align  itself  more  axially 
than  for  a  shorter  length  duct,  and  therefore,  the  momentum  drag  is  greater. 

In  order  to  overcome  the  larger  total  drag,  the  thrust  vector  of  the 
vehicle  must  be  tilted  sufficiently  co  that  the  component  of  thrust  in  the 
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horizontal  direction  exactly  cancels  the  total  drag  in  steady-state  flight. 
Although  the  short  chord  duct  yields  the  more  desirable  smaller  tilt  angles 
other  characteristics  such  as  the  performance  of  the  long  chord  duct 
would  hlso  be  considered  in  any  design  (Ref.  2). 


Static.  Stability 

1 .  The  E  Coefficient 

As  mentioned  previously  in  the  DEVELOPMENT,  the  fundamental 
mathematical  condition  for  static  stability  is  that  the  E  coefficient  (the 
term  independent  of  X  in  the  stability  quartic)  be  positive. 

The  E  coefficient: 


/  — 


^  ~  C*  s) 


s4  oS" 


(38) 


is  a  function  of  the  two  static  stability  partial  derivatives  m  and 

C  ii 


m 


a 


The  term  m  is  referred  to  as  the  static  velocity  stability  and 

C  ^ 

m  the  static  angle  of  attack  stability. 


O' 


An  aircraft  is  statically  stable  with  velocity, if,  when  acted  upon  by 
a  pure  velocity  disturbance,  it  creates  a  moment  about  the  center  of  gravity 
tending  to  restore  the  aircraft  to  the  initial  velocity.  Similarly,  an  air¬ 
craft  is  statically  stable  with  angle  of  attack,  if,  when  acted  upon  by  a 
pure  angle  of  attack  disturbance,  it  creates  a  moment  about  the  center  of 
gravity  tending  to  restore  the  aircraft  to  the  initial  angle  of  attack.  The 

Q 

criterion  for  static  velocity  stability  requires  that  m  be  greater  than  zero 

C  ^ 

while  static  angle  of  attack  stability  requires  that  m  be  less  than  zero. 

O' 
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Generally  the  Flying  Platform  has  too  much  velocity  stability  in 
hovering  and  forward  flight.  However  the  velocity  stability  does 
decrease  as  flight  speed  increases  and,  depending  on  c.g.  iocation,  may 
even  become  unstable  by  changing  sign. 

Q 

Wi  h  respect  to  the  angle  of  attack  stability,  in  hovering  m 

a 

is  equal  to  zero  and  the  Flying  Platform  is  neutrally  stable  with  angle  of 

c  c 

attack.  In  forward  flight,  m  has  the  wrong  sign,  that  is,  m  increases 

a  a 

positively  as  flight  speed  increases  so  that  the  platform  is  statically  unstable 
with  angle  of  attack.  Neither  of  these  criteria,  taken  by  itself,  is 
adequate  to  predict  the  static  stability  of  the  Flying  Platform.  The  commonly 

Q 

used  static  stability  criterion  of  m  <  0  for  airplanes,  is  not  sufficient  for 

C  a 

this  case  since  it  assumes  m  is  equal  to  zero. 

For  an  aircraft  operating  in  a  flight  regime  or  configuration  where 

0 

m  is  not  equal  to  zero  }  the  general  criterion  for  static  stability  must 
^  G  G 

include  the  effects  of  both  m  and  m  as  exemplified  by  the  E  coefficient. 

§1  a 

Static  instability  results  in  a  dynamic  divergence  whenever  this  coefficient 
becomes  negative.  .  In  the  speed  regime  considered  for  this  Flying  Platform, 

c  c 

all  the  terms  in  the  E  coefficient  are  always  positive  except  m^  and  T  . 

For  the  long  chord  duct,  the  term  ^m^  is  generally  positive  and  large, 

but  can  be  reduced  and  even  made  negative  in  certain  forward  flight  speeds 

0 

by  raising  the  center  of  gravity.  The  term  T  becomes  negative  rather  quickly 

0  M  0  0 

and  since  m^  is  always  positive  the  quantity  -  T^  m^  is  usually  a 

positive  increment  to  the  E  coefficient.  From  this  it  is  evident  that  the  E 

coefficient  of  the  long  chord  duct,  in  the  region  of  speeds  considered,  will 

probably_not  change  sign  unless  the  center  of  gravity  is  raised  sufficiently 

Q 

so  that  m  can  become  a  negative  value. 

M 
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Also,  It  is  interesting  to  note  for  this  case  that  the  larger  the  posl- 

c  c  c 

tive  (unstable)  value  of  m  the  laraer  the  value  of  tne  quantity  -  T  . 

a  1  fi  a 

This  becomes  a  larger  positive  increment  to  the  E  coefficient  and  therefore 

c 

enhances  the  static  stability.  However,  large  positive  m  has  a  simul- 

a 

taneous  adverse  effect  on  the  other  coefficients  of  the  stability  quartic. 

C 

For  example,  for  large  m  ,  positive  static  stability  may  exist  while  in 

a 

all  probability  the  C  coefficient  will  have  already  changed  sign,  indicating 
a  serious  instability.  See  Table  I. 

c 

For  the  short  chord  duct,  the  values  of  m  are  considerably 
smaller  and  become  negative  more  quickly  than  the  long  chord  duct  as 

Q 

forward  flight  speed  is  increased.  Also,  remains  positive  for  the 
speeds  considered.  The  first  term  in  the  E  coefficient  becomes  a  negative 

Q 

increment  as  soon  as  m  becomes  negative.  The  second  term  is  always 

^  C 

positive  so  that  E  changes  sign  as  sopn  as  the  value  of  m  has  reduced 

c  c  c  c  ^  c 

to  a  small  positive  number  where  T  m  >  m  T  .  Whenever  m 

li  a  li  u  ju 

becomes  negative,  E  has  already  become  negative  for  this  case. 

Although  it  is  desirable  to  have  a  positive  static  stability,  insta.- 
bilities  due  to  E  changing  sign  are  associated  with  velocity  changes  and 
are  probably  not  as  serious  as  the  instabilities  associated  with  angle  of 
attack  changes  which  are  indicated  by  other  coefficients  of  the  stability 
quartid  changing  sign. 

2.  The  Moment  to  Trim  Versus  Velocity  Diagram 

The  static-trim  curve,  in  which  the  moment  produced  by  the 
machine  in  trim  equilibrium  is  plotted  versus  forward  speed,  can  also 
be  used  to  investigate  the  static  stability  of  the  Flying  Platform. 

In  Fig.  19,  the  M-V  curves  are  plotted  for  both  ducts  for  the 
various  center  of  gravity  locations  by  using  the  moment  equation: 


Lily  .=  Mc/  +T-  xt  ■  H  ~  ° 


c 


(4  c) 
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where  T  •  x^,  is  the  control  input  of  the  pilot,  and  the  moment  produced 
by  the  vehicle  including  the  drag  of  the  pilot  is: 


(61) 


Inspection  of  the  curves  of  this  moment  plotted  versus  velocity 
indicates  that  the  value  of  the  moment  required  to  trim  can  be  reduced 
by  raising  the  center  of  gravity  position. 

The  static  stability,  or  slope  of  the  M-V  trim  curve,  is  the  total 

.  In  all  cases  presented,  the  slope  starts  as  positive, 


derivative  <&L 
Oi  v 

reaches  a  zero  value,  and  then  becomes  negative.  The  criterion  for 
static  stability  is 


dV 


>  0. 


Also  plotted  on  the  M-V  diagram  is  the  horizontal  line  drawn 
at  200  foot-pounds  representing  the  approximate  maximum  control  moment 
attainable  by  the  pilot  for  trim.  Presuming  that  a  typical  200-pound  pilot, 
with  his  c.g.  located  at  3^  feet,  can  lean  forward  20  degrees  relative  to 
the  platform,  the  maximum  control  moment  he  can  apply  would  be  approxi¬ 
mately: 

/i  =  (2.  00)( 3.  3  3)(s<*  20  °) 

/VJ/. 


MP  =>  200  ft.'  lira. 


(62) 


By  convention,  the  distance  x  is  taken  as  minus  when  the  c.g. 
is  moved  forward  in  forward  flight.  The  scale  for  plotted  opposite  M„ 
then,  in  effect,  represents  control  position  for  trim. 

In  forward  flight,  the  maximum  forward  displacement  is  equal 
to  approximately  -0.36  feet  for  200  foot-pounds.  The  intersections  of 
the  curves  and  the  horizontal  line,  (denoted  by  circles),  represent  the 


V 


51 


maximum  forward  speed  for  which  the  pilot  is  able  to  trim  the  moment 
produced  by  the  Platform. 


3 


The  short  vertical  dashes  on  each  curve  represent  the  points  where 


sr 


equals  zero  and  beyond  which 


becomes  negative .  These 


V  .  ,  dv 

points  are  exactly  the  sppeds  at  which  the  E  coefficient  changes  sign  and 
becomes  negative.  Actually  the  speeds  beyond  the  vertical  dash  location 
(»•<>).  represent  flight  in  a  speed  regime  where  the  Platform  is 
statically  unstable. 

For  a  statically  stable  aircraft,  the  control  position  for  trim  always 
moves  forward  (piiot  must  be  positioned  further  forward  to  trim  the  Platform) 
as  higher  trim  flight  speeds  are  achieved. 

Notice  on  the  diagram  that  the  control  position  for  trim  ,  starts 
to  reduce  beyond  these  vertical  dash  points.  In  this  region,  the  control 
position  Xj,  for  trim  (pilot  position)  is  located  at  successively  smaller 
negative  values  for  greater  trim  speeds  after  having  reached  the  maximum 
negative  value  at  equal  to  zero.  This  should  not  be  thought  of 

as  a  region  of  so  called  "control  reversal",  in  the  normal  sense  of  the 
words,  even  though  the  control  moves  aft  to  trim  each  successively  higher 
flight  speed.  Usually  the  term  "control  reversal"  implies  quite  a  different 
meaning  with  different  effects  and  is  most  frequently  associated  with  the 
reversed  control  effectiveness  due  to  aeroelastic  phenomena.,  when  the 


control  derivative , 


e$C 

-v rvrnk.  '  changes  sign  (6  is  the  control  deflection). 

d&c  c 


With  this  meaning,  for  example,  whenever  a  pilot  attempts  to  pull  the  nose 
of  an  airplane  up  in  order  to  reduce  speed,  (control  moved  aft),  the  aircraft 
actually  noses  down  and  tends  to  gain  speed.  Fortunately,  this  is  not  what 
occurs  for  the  Flying  Platform  if  it  is  operated  in  a  statically  unstable 
configuration,  or,  for  that  matter,  ♦■o  any  airplane  which  happens  to  be 
operating  in  a  statically  unstable  configuration.  In  order  to  go  faster  or 
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move  the  nose  down,  It  Is  still  initially  necessary  to  move  the  control 
forward  (lean  forward)  to  obtain  the  desired  increase  in  speed.  What  is 
different  Is  that  once  the  desired  higher  speed  is  achieved,  and  the  vehi¬ 
cle  is  again  in  trim,  the  pilot  will  find  the  control  positioned  forward  less 
than  it  was  for  the  slower  trim  position  (pilot  leaning  forward  less  than 
previous  trim  position).  In  order  to  accelerate  an  aircraft  in  a  given 
direction,  the  control  will  always  be  initially  moved  in  the  same  sense 
regardless  of  whether  the  aircraft  is  statically  stable  or  unstable. 

To  slow  down,  the  pilot  initially  leans  back,  but  when  the  aircraft 
is  again  trimmed  at  the  desired  slower  flight  speed ,  in  the  region  where 

OV  iS  ne9atlve'  he  wm  And  he  is  leaning  further  forward  than  before. 
As  seen  by  the  curves  in  the  M-V  diagram,  raising  the  c.g.  does 
reduce  the  control  moment  required  for  trim  but,  unfortunately,  also  moves 

this  point  (  E  =  -ffi-  --  0)  to  successively  lower  forward  speeds,  thereby 
reducing  the  statically  stable  flight  envelope. 

cLM  ThS  aSt6riSks  °n  the  dlagram  indicate  the  maximum  points  where 
d  /  13  equai  t0  zer0  -  for  eacil  duct,  commensurate  with  the  200 

foot-pound  control  moment.  The  difference  between  the  two  radically 

different  duct  designs  appears  to  be  on  the  order  of  only  six  feet  per 
second: 

Short  Chord  Duct  »  38  feet  per  second 

Long  Chord  Duct  vmax  ^  32  feet  per  second 

Other  flying  qualities  of  a  machine  similar  in  nature  to  the 
subject  vehicle  are  discussed  in  Ref.  3. 

Dynamic  Stability 

The  fundamental  requirement  for  stable  dynamic  response  is  that 

all  the  coefficients  of  the  stability  quartic  and  Routt, 's  discriminant  be 
positive. 


M- 


W ,y„ 
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If  ail  the  coefficients  are  initially  positive,  the  first  indication  of 
an  instaoility  will  normally  be  a  divergence  when  E  changes  sign,  or  an 
unstable  oscillation  when  Routh's  discriminant  becomes  negative. 

The  roots  of  the  characteristic  equations  evaluated  for  the  original 
c.g.  location  of  each  duct  at  equal  to  0.05  and  0.10  are  shown  in  Tables 
I  and  II. 

It  should  be  remembered  that  the  following  analysis  investigates 
the  control-fixed  or  open-loop  response  of  the  Flying  Platform  configurations. 

The  open-loop  response  is  the  resulting  mode  of  motion  of  the 
Flying  Platform  when  the  pilot  remains  rigid  and  applies  no  control  inputs 
throughout  the  transient  response.  The  closed-loop  response  is  the 
resulting  mode  of  motion  of  the  Platform  when  the  pilot  acts  as  a  feedback 
by  providing  controlling  inputs  to  the  transient  response.  Although  the 
Flying  Platform  is  unstable  in  the  open-loop  it  will  be  stable  closed- 
loop  if  the  pilot  can  correctly  control  the  machine  as  he  desires. 

The  open-loop  response  of  the  Flying  Platform  is  investigated  to 
determine  the  effects  on  the  dynamic  stability  of  vertical  center  of  gravity 
location,  duct  design  changes  and  artificial  stabilization. 

In  this  Report  this  may  be  thoroughly  accomplished  by  referring 
to  the  various  figures  of  the  stability  diagrams,  boundaries  of  the  quartic 
and  the  associated  root  locus  plots. 

On  the  figures  labeled  "Boundaries  of  the  Stability  Quartic"  and 

c  c 

"Stability  Diagram",  the  symbols  represent  the  values  of  m  and  m 

M  Q 

for  each  particular  center  of  gravity  location.  The  dashed  line  connecting 
these  circles  represents  the  movement  of  these  points  on  the  diagram, 
as  the  center  of  gravity  location  is  varied. 


4*.  , 
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For  the  different  cases: 

1 .  Long  and  Short  Chord  Duct  Hovering  p  =  V  =  Q 

The  dynamic  analysis  of  the  hovering  case  of  the  subject  Flying 

Platform  can  best  be  obtained  by  referring  to  the  root  locus  diagram  Fig.  7 
a  and  b. 


The  case  is  somewhat  simplified  by  the  fact  that  we  deal  with  a 
cubic  equation  and  is  equal  to  zero. 


For  the  hovering  long  chord  duct  (Fig.  7b),  when  gm  Is  equal  to 
zero,  the  three  poles  represent  two  convergences  (point  80  and  point  8J) 
and  a  zero  root  (point  82) /  As  the  gain  of  Cm  is  Increased  positively, 

the  pole  at  point  80  moves  toward  point  86  and  continues  as  a  convergence 
for  all  values  of  positive  m 
C  M 

Also,  as  m^  is  Increased  positively,  the  two  poles  at  point  81 
and  point  82  move  toward  each  other  on  the  negative  real  axis  until  they 
become  equal  roots  at  the  break  point  fcolnt  83).  Further  increases  in  dm 
cause  the  equal  roots  to  break  away  from  the  negative  real  axis  and  become 
a  complex  pair,  indicating  a  stable  oscillation  of  Increasing  frequency 
and  decreasing  stability.  At  point  84,  Cm  has  Increased  to  the  value  at 
which  the  oscillation  becomes  neutrally  damped.  Increasing  further 


makes  the  oscillation  unstable  (point  85)_.  The  oscillation  continues  to 
Increase  In  instability  and  frequency  as  ^  is  varied  to  positive  Infinity 

The  short  chord  duct  In  hovering  (Fig.  7a)  has  a  similar  analysis 
and  is  numbered  accordingly,  the  only  difference  being  the  respective 
gains  for  the  particular  modes  ot  motion. 

In  order  to  complete  the  set  of  stability  diagrams.  Fig.  8  is 
presented  showing  the  various  gains  of  dm  for  the  different  center  of 
gravity  locations  together  with  the  various  modes  of  motion.  This 

diagram  is  considerably  different  from  the  other  stability  diagrams  since 
only  the  values  of  m  need  to  be  shown. 
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For  the  long  chord  duct  m  changes  approximately  0.0028  per 

M  Q 

inch  of  c.g.  movement.  Stable  oscillations  start  when  m  is  equal  to 

C  ^ 

0.0003,  which  is  the  value  of  m  at  the  root  locus  break  point  (point  83). 

M  c 

Unstable  oscillations  start  when  m  is  equal  to  0.0056,  which  is  the 
C  ^ 

value  of  m  at  point  84  on  the  root  locus  diagram.  All  the  center  of 

M 

gravity  locations  are  well  inside  the  unstable  region. 

Q 

For  the  short  chord  duct,  m  changes  approximately  0.001  per 

/i  C 

inch  of  c.g.  movement  Stable  oscillations  start  when  m  is  equal  to 

C  M 

0.000C75,  which  is  the  value  of  m  at  the  root  locus  break  point  (point  93). 

C  ^ 

Unstable  oscillations  start  when  m  is  equal  to  0.0024,  the  value  of 

C  M 

m^at  point  94  on  the  root  locus  diagram.  Again,  all  the  c  g.  locations 

are  well  inside  the  unstable  region. 

For  the  specific  case  of  the  long  chord  duct  in  hovering,  the  roots 
of  the  characteristic  equation  (Table  I)  for  the  c.g.  in  the  original  location 
indicate  a  convergent  mode  and  an  unstable  oscillation.  The  period  of 
this  oscillation  is  approximately  four  seconds  with  a  time  to  double 
amplitude  of  about  one  second. 

For  the  short  chord  duct  in  hovering,  the  roots  of  the  characteristic 
equation  (Table  II)  for  the  c.g.  in  the  original  location  indicate  a  conver¬ 
gent  mode  and  an  unstable  oscillation.  The  period  of  this  oscillation  is 
approximately  six  seconds  with  a  time  to  double  amplitude  of  2.23  seconds. 

Both  ducts  are  rather  unstable  with  regard  to  the  open-loop  response, 
however,  the  short  chord  duct  is  noticeably  the  better  of  the  two. 

As  mentioned  in  the  DEVELOPMENT,  it  is  interesting  to  investigate 
the  effects  of  artificial  stabilization  on  the  transient  response  of  these 
vehicles . 

The  effect  of  such  devices  on  the  behavior  of  the  system  is  shown 
by  root  locus  plots  for  several  different  cases. 
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2  .  Long  and  Short  Chord  Duct  Artificial  Stabilization  Added 

Hovering  [X  =  0 

For  the  short  chord  duct  (c.g.  in  the  original  location)  when  K 
varies  from  zero  to  positive  infinity  reference  is  made  to  Fig.  9a. 

The  artificial  stabilization  adds  two  zeros  to  the  hovering  root  locus 
diagram,  one  at  the  origin  and  the  other  close  by  in  the  left,  half  plane. 
Arbitrarily  selecting  T  equal  to  3.333,  an  additional  pole  is  located  at  -0.3. 

When  K  is  equal  to  zero  ,  the  poles  located  at  point  100  indicate  an 
unstable,  lightly  damped  oscillation,  while  the  poles  at  point  101  and  point 
102  represent  two  convergences.  As  K  is  increased,  the  poles  (points  101  and 
102)  come  together  on  the  real  axis  at  the  break  point  103  and  then  bre^k 
away  from  the  real  axis  to  become  a  stable  oscillation  of  increasing 
frequency.  Meanwhile  the  unstable  oscillation  becomes  stable  with  no 
appreciable  change  in  frequency  (point  104).  However,  when  K  is  increased 
further,  the  oscillation  remains  stable  but  decreases  in  frequency  (point  106) 
until  the  complex  pair  become  real  roots  which  have  limiting  values  at  the 
zeros  on  the  real  axis  (point  106). 

If  we  make  T  larger,  for  example  equal  to  50,  so  that  the  pole  is 
located  between  the  zeros  on  the  real  axis,  the  root  locus  appears  as  in 
Fig.  9b. 

When  K  is  equal  to  zero,  the  poles  (point  112  and  point  109)  rep¬ 
resent  convergences  and  the  complex  pair  (point  110)  an  unstable  oscillation. 
As  K  is  increased,  the  poles  representing  the  convergences  (point  112  and 
point  109)  move  toward  and  terminate  at  the  zeros  (point.  108  and  point  107). 
The  pole  at  point  109  becomes  a  zero  root  when  K  reaches  infinity,  whereas 
the  other  pole  (112)  always  represents  a  convergence. 

* 

The  unstable  oscillation  (point  110)  increases  in  frequency  and 
becomes  neutrally  stable  at  point  111.  As  K  is  increased,  the  oscillation 
remains  stable  but  the  frequency  continues  to  increase. 

As  suspected,  depending  on  the  value  of  the  time  lag,  the  addition 
of  a  lagged  rate  feedback  device  is  able  to  stabilize  the  hovering  platform 
to  various  degrees . 
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For  pure  rate  feedback,  when  T  in  is  selected  equal  to  zero  (T  = 

g  =  -  =  0),  the  extra  pole  is  located  at  negative  infinity  (a  =  -oo). 

As  K  is  increased  positively,  the  poles  representing  the  unstable  oscillation 

(point  100)  move  directly  to  the  break  point  i06  in  the  left  half  plane  instead 

of  looping  around  (points  104  and  105)  as  in  Fig.  9a.  For  this  case  the  pole 

representing  the  convergence  at  point  102  stays  in  the  left  half  plane  and 

moves  to  the  left  as  K  is  increased. 

For  pure  attitude  feedback,  when  T  is  selected  as  infinitely  large 

(T  =  ;  =  k  *  oo),  the  extra  pole  (point  109)  is  located  at  the  origin  (a  =  0) 
a  u 

and  cancels  the  zero  (point  107)  located  there.  The  root  locus  for  this  case 
appears  very  similar  to  Fig.  9b. 

Thus,  the  hovering  Flying  Platform  may  also  be  stabilized  by  either  pure 
rate  or  attitude  feedback  or  a  combination  of  both . 

The  analysis  of  adding  artificial  stabilization  to  the  hovering  long 
chord  duct  is  qualitatively  the  same  as  that  of  the  short  chord  duct. 

3.  Long  Chord  Duct  ii  =  0.05  V  =  16  mph. 

Referring  to  the  boundaries  of  the  quartic  (Fig.  10A);  it  is  seen 

Q 

that  center  of  gravity  movement  has  a  primary  effect  on  the  value  of  m 

C  ** 

but  only  changes  m  slightly.  Actually,  on  this  diagram,  each  inch 

O'  ~ 

C 

of  c.g.  movement  accounts  for  approximately  0.0031  change  in  m^  and 

0.000014  change  in  .  It  is  evident  that  all  c.g.  locations  are 

a 

on  the  right  side  of  the  C  =  0  line  and  therefore  the  C  coefficient  is  always 
negative.  If  the  c.g.  is  raised  more  than  six  inches  above  the  original 
position,  the  D  coefficient  changes  sign  (point  1),  and  if  raised  by  12 

inches  (point  2),  the  E  coefficient  changes  sign. 

* 

Referring  to  the  stability  diagram  (Fig.  1 0 B) ;  It  is  noted  that  most 
of  the  c.g.  locations  fall  within  the  region  of  unstable  oscillations.  If 
the  c.g.  is  raised  about  nine  inches  (point  3)  the  unstable  oscillations 
cease  and  divergences  occur.  Continuing  to  raise  the  c.g.  (point  4)  causes 

9 

E  to  change  sign  so  that  even  the  static  stability  is  now  negative. 
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A  better  appreciation  of  the  modes  of  motion  can  be  obtained  from 

c 

the  associated  root  locus  plots.  Note  that  the  value  of  m  on  the 

a 

stability  diagram  is  approximately  equal  to  0.002  and  corresponds  approxi¬ 
mately  to  the  value  investigated  in  Fig.  10-I-c. 

C 

When  is  equal  to  zero,  the  poles  are  located  at  point  5  and 
indicate  the  presence  of  two  convergent  and  two  divergent  modes.  The 

Q 

value  or  gain  of  m  Increases  from  zero,  at  the  poles,  to  infinity  along 

Q 

the  direction  of  the  arrows.  As  m  is  increased  to  the  value  at  point  3 

M 

(the  break  point  on  the  root  locus  and  corresponding  to  point  3  in  Fig.  10B), 
an  unstable  oscillation  starts  while  the  two  other  roots  move  along  the 

arrows  and  remain  convergent  modes.  Continuing  to  increase  the  gain  of 

C 

m^  to  larger  positive  values  along  point  6,  the  oscillation  increases 
in  frequency  and  instability.  The  two  other  poles  remain  convergences 

Q 

regardless  of  the  positive  gain  of  m  . 

C 

The  portion  of  the  c.g.  travel  at  which  m  is  negative  may  be 

*  ^  C 

investigated  by  referring  to  Fig.  10-IV-c.  Again  at  zero  m^  ,  there  are 
the  aperiodic  modes,  two  convergences  and  two  divergences  (point  5).  As 

C 

the  gain  of  is  increased  negatively,  it  is  seen  that  a  stable  oscillation 
starts  (point  8)  and  one  divergent  mode  becomes  convergent.  However, 
the  predominating  mode  is  the  ever-increasing  divergence  (point  7),  as 

c  c 

m  is  increased  negatively  so  that  the  overall  mode  of  motion  for 
increasing  negatively  is  an  ever-increasing  divergence. 

The  open-loop  response  in  this  case  is  a  serious  instability.  If 
the  roots  calculated  for  the  case  with  the  original  c.g,  location  are 
considered  (Table  I),  it  is  seen  that  at  ii  =  .  05  and  V  =  16  mph,  there  is 
an  unstable  oscillation  with  a  period  of  8.57  seconds  and  a  time  to 
double  amplitude  of  only  a  little  over  two  seconds.  In  other  words,  the 
oscillation  would  double  amplitude  in  about  one -quarter  of  the  cycle. 
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4.  Long  Chord  Duct  /i  =  0.10  V  =  34  mph 

Referring  to  the  boundaries  of  the  quartic  (Fig.  11A);  it  is  seen 

C 

that  m  changes  approximately  0.0042  per  inch  of  c.g.  movement  while 

C  ^ 

m  changes  0,00014  per  inch. 
a 

Starting  with  the  center  of  gravity  lowered  six  inches  from  the  ori¬ 
ginal  location,  it  is  noted  that  the  c.g.  is  on  the  positive  side  (point  11) 

of  the  C  =  0  line.  As  the  c.g.  is  raised  toward  the  original  location,  the 

C  C 

C  coefficient  changes  sign  (point  12)  and  m  decreases  while  m 

p.  a 

increases.  Continuing  to  raise  the  c.g.  further,  the  D  coefficient  changes 
sign  (point  13).  When  the  center  of  gravity  is  raised  six  inches  (point  14), 

c 

the  line  RLBP  has  been  crossed  and  m  is  negative.  If  the  c.g.  is  raised  a 

P 

few  more  inches  (point  15)  the  E  coefficient  changes  sign. 

Referring  to  the  stability  diagram  (Fig.  1  IB); the  points  again 

are  located  principally  in  the  region  of  unstable  oscillations.  The  unstable 

oscillations  cease  if  the  c.g.  is  raised  approximately  six  inches  to  point 

14  at  which  location  the  RLBP  line  has  been  crossed  and  the  dynamic 

instability  takes  the  form  of  a  divergence.  If  the  c.g.  is  raised  a  few 
C 

more  inches,  m  continues  to  increase  negatively  and  at  point  15,  E 

P 

changes  sign^so  that  now  the  configuration  is  even  statically  unstable. 

For  a  closer  look  at  the  characteristic  modes  of  motion  reference 

is  again  made  to  the  appropriate  root  locus  plots.  From  the  stability 

C 

diagram  note  that  m  ranges  in  value  from  approximately  0.001  to  0.004. 

a 

The  root  locus  plots  for  both  these  values  are  shown  in  Figs.  11-I-d 
and  11-I-e. 

Fig.  11-I-d  represents  the  locus  when  moving  vertically  up  the 

c  c 

line  m  equal  to  0.001  in  the  stability  diagram  (Fig.  1  IB)  when  m 
«  P 

varies  from  zero  to  positive  infinity.  The  initial  position  of  the  poles, 

C  C 

when  m  equals  0.001  and  m  equals  zero,  indicates  that  two  coh¬ 
o'  p, 

vergences  (point  16)  and  an  unstable  oscillation  (point.  17)  exist.  As 


! 
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the  gain  of  m  is  increased,  (equivalent  to  moving  upward  along 
C  P 

m  equal  to  0.001  in  the  stability  diagram)  by  moving  in  the  direction 
a 

of  the  arrows  in  the  root  locus  diagram,  two  poles  always  remain  con¬ 
vergences,  one  with  ever  increasing  value  (point  18)  and  the  other 
with  decreasing  value  and  a  limit  at  the  zero  (point  19).  The  other  two 
poles  (point  17)  continue  as  unstable  oscillations  with  ever-increasing 

Q 

frequency  a qd  instability  (point  20)  as  m  is  increased. 

C  ^ 

For  m  equal  to  0.004  (Fig.  11-I-e),  the  picture  is  much  the  same, 
®  C 

except  that  for  m  equal  to  zero,  the  four  poles  represent  two  conver- 

P  Q 

gences  (point  21)  and  two  divergences  (point  22).  As  the  gain  of  m^ 
is  increased  the  convergences  behave  as  before.  However,  the  two  poles 
representing  the  divergences  (point  22)  move  toward  each  other  and  become 
equal  roots  at  point  23,  the  root  locus  break  point.  This  point  (point  23) 

Q 

is  also  shown  on  the  stability  diagram.  As  the  gain  of  m^  is  increased 
further,  the  two  roots  split  as  shown  (point  24)  and  become  unstable 
oscillations  of  increasing  frequency  and  instability. 

In  order  to  cover  the  remaining  center  of  gravity  locations ,  (point 

c 

14  and  point  15  in  Fig.  11B),  the  root  locus  plots  for  m  equal  to  the  same 

C  »  ” 

values,  but  varying  from  zero  to  negative  infinity,  are  analyzed 

(Fig.  11-IV-d  and  Fig.  11-IV-e). 

Q 

In  Fig.  11-IV-d,  the  four  poles  at  equal  to  zero  represent 
two  convergences  (point  25)  and  an  unstable  oscillation  (point  26).  As 

c 

the  gain  of  is  increased  negatively,  the  two  convergent  modes 
(point  25)  move  toward  each  other  in  the  left  half  plane  to  the  break  point 

Q 

(point  27).  As  m  is  increased  more  negatively,  they  break  away  from  the 

P 

real  axis  and  become  a  complex  pair  representing  a  stable  oscillation 
of  increasing  frequency  and  decreasing  damping  ratio.  In  the  meantime, 

c 

as  m  is  being  increased  from  zero  to  negative  infinity,  the  unstable 

P 


'X4& 
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oscillation  (point  26)  is  reducing  frequency  to  the  point  where  the  complex 
pair  of  roots  (point  26)  become  divergences,  or  equal  real  roots  by  joining 
on  the  real  axis,  (point  28).  One  divergent  pole  formerly  at  point  28 
moves  across  the  origin  and  becomes  a  convergence  with  a  limiting  value 
at  the  zero  (point  29)  in  the  left  half  plane.  The  other  divergent  pole, 
formerly  at  point  28,  moves  to  the  right,  continuing  as  a  divergence  of 
ever-increa3_ing  instability  (point  30). 

Q 

For  m  equal  to  0.004  (Fig.  11-IV-e),  the  poles  behave  in  much  the 
a 

same  fashion  with  the  exception  that  the  unstable  oscillation  at  point  26 
becomes  two  unequal,  positive,  real  roots,  representing  divergences  when 

c 

m  equals  zero.  These  poles  move  in  much  the  same  way  as  the  posi- 

M  Q 

tive  real  roots  (point  2  8)  move  as  m  increases  negatively.  Note  that  in 

M 

both  cases  there  is  always  an  instability  present,  either_in  the  form  of 

Q 

a  divergence  or  an  unstable  oscillation  for  the  range  of  m  from  zero 
to  negative  infinity. 

Again,  the  open-loop  response  for  this  case  is  a  serious  insta¬ 
bility.  Considering  the  roots  calculated  for  the  case  of  the  original  c.g. 
location  (Table  III),  for  /i  =  0.10  and  V  =  34  mph,  two  convergences  and 
an  unstable  oscillation  exist.  The  predominating  mode  of  motion  for  this 
case  is  the  unstable  oscillation  with  a  ten  second  period  and  a  two  second 
time  to  double  amplitude.  Therefore,  the  aircraft  doubles  amplitude  in 
approximately  1/5  of  a  cycle. 


5,.  Long  Chord  Duc  t  Artificial  Stabilization  A^ded  p.  =  0.10  V=  34mph 
Arbitrarily  choosing  to  investigate  the  effects  of  artificial  stabi¬ 
lization  on  the  dynamic  stability  associated  with  the  vertical  center  of 
gravity  located  at  the  original  location,  an  instability  exists  in  the  form 
of  an  unstable  oscillation  (Fig.  12,  points  140  and  150).  Taking  the  case 
of  pure  attitude  feedback  (Fig.  12a);  one  pole  cancels  with  a  zero  at  the 
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origin  (T  equals  infinity),  the  remaining  poles  represent  two  convergences 
(point  141  and  point  142)  and  the  unstable  oscillation  (point  140).  As  K 
is  increased,  the  unstable  oscillation  (point  140)  becomes  neutral  (point 
146)  and  finally  stable,  with  the  oscillation  reaching  a  limiting  value  at 
the  zeros  (point  143).  The  two  poles  at  point  141  and  point  142  join  at 
point  144  and  break  away  from  the  negative  real  axis  to  become  a  stable 
oscillation  of  increasing  frequency  (point  145)  as  K  is  increased. 

For  the  case  of  pure  rate  feedback  (Fig.  12b),  a. zero  is  located  at 
the  origin  (T  equals  zero),  and  the  other  poles  and  zeros  are  located  in 
the  same  position  as  in  Fig.  12a.  The  unstable  oscillation  (point  150) 
and  the  two  convergent  poles  (points  151  and  152)  still  exist  as  before. 

However,  as  the  gain  K  is  increased,  the  poles  of  the  unstable 
oscillation  (point  150)  move  into  the  left  half  plane  (point  156)  and  be¬ 
come  a  stable  oscillation  with  a  limiting  value  n t  the  zeros  (point  154). 
Notice  that  they  move  in  a  more  direct  line  contrasted  to  the  former  case 
where  the  loci  loop  around  (point  146)  when  moving  toward  the  zero  (143). 

To  investigate  intermediate  values  of  T  (pure  rate  feedback  with 
various  degrees  of  lag),  refer  to  Fig.  12a.  Locating  a  zero  at  the  origin 
w<th  a  pole  very  close  to  it  but  in  the  left  half  plane,  the  loci  are 
essentially  unchanged.  As  this  pole  moves  to  the  left  or  away  from  origin 
(T  decreasing),  the  loci  between  point  14  0  and  point  143  tend  to  describe 
paths  that  become  more  and  more  direct  until  they  appear  exactly  like 
the  loci  between  point  150  and  154  when  the  moving  pole  reaches  negative 
infinity.  The  other  loci  for  point  141  and  142  remain  similar  except  that 
the  break  point  (point  144)  continues  to  move  to  the  left  as  the  moving 
pole  passes  the  position  point  142.  The  pole  at  point  142,  in  this  case, 
always  moves  to  the  right  and  terminates  at  the  zero  located  at  the  origin. 

From  this  analysis,  it  is  evident  that  a  pure  rate  feedback  with 
various  degrees  of  lag  can  be  used  to  stabilize  the  system  and  exert  some 
control  on  the  frequency  and  damping  ratio  of  this  particular  configuration. 


t 
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6.  Short  Chord  Duct  /i  =  0. 05  V  =  18  mph 

Referring  to  the  boundaries  of  the  quartic  (Fig.  13A); 


m 


changes  approximately  0.002  per  inch  of  vertical  c.g.  movement  while 

<5 

m  changes  very  slightly,  say  0.00001  per  inch  movement.  All  center 
®  C 

of  gravity  locations  have  positive  m  and  are  located  to  the  left  (positive 

a 

side)  of  the  C  =  0  line  so  that  the  C  coefficient  in  the  stability  quartic 
is  alv/ays  positive  for  this  case.  However,  if  the  c.g.  is  raised  almost 
six  inches  above  the  original  location  (point  34),  the  E  coefficient  changes 
sign  so  that  a  dynamic  instability  appears  in  the  form  of  a  divergence, 
due  to  the  static  instability.  Raising  the  c.g.  slightly  further  (point 

Q 

35),  m^  becomes  negative  and  an  additional  increment  of  upward  travel 
causes  the  D  coefficient  to  change  sign  as  well  (point  35). 

Referring  to  the  stability  diagram  (Fig.  13B);  the  center  of  gravity 
located  six  inches  below  the  original  location  lies  in  a  region  of  unstable 
oscillations  (point  31).  As  the  center  of  gravity  is  raised  past  the  original 

Q 

c.g.  location  (point  32)  (  m  decreasing)  it  crosses  the  line  R  =  0  (point 

33),  thereby  ending  the  unstable  oscillation, and  passes  through  a  small 

portion  of  the  stable  region  before  crossing  the  E  =  0  line  (point  34) 

C 

into  the  divergent  region.  If  the  c.g.  is  raised  further,  changes 
sign  (point  35) . 

O 

To  detefmine  the  characteristic  modes  of  motion  for  this  case, 
reference  is  again  made  to  the  appropriate  root  locus  plot.  From  the  stability 

Q 

diagram,  m  is  in  the  range  of  approximately  0.0002.  The  closest 

a  Q 

value  in  the  root  locus  plots  is  m  equal  to  0.0005.  Analyzing  the 

a 

plot  (Fiq.  13-I-b),  which  closely  resembles  the  character  of  the 

Q 

motion  for  the  exact  value  of  m  ,  the  following  statements  can  be  made: 

C  ^  C 

When  m  equals  0.0005  and  m  equals  zero,  one  pole  (point  36) 

O'  ii 

indicates  a  convergence,  the  second  (point  37)  a  divergence  and  the 
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remaining  two  (point  38)  a  stable  oscillation.  This  is  approximately 
the  situation  existing  at  point  35  on  the  stability  diagram  (Fig.  13B) 

Q 

for  the  c.g.  location  six  inches  above  the  original.  As  m  is  increased 

M 

from  zero  to  positive  infinity,  the  pole  at  point  36  continues  as  an  increas¬ 
ing  convergence  (point  48),  i.e.  ,  this  transient  mode  tends  to  die  out 
faster.  The  divergent  pole  at  point  37  moves  from  the  right  half  plane, 
across  the  origin  (point  34),  becoming  a  convergence  in  the  left  half  plane 

Q 

and  reaching  a  limiting  value  at  the  zero  (point  49)  when  m  is  increased 

M 

to  positive  infinity. 

With  respect  to  the  two  poles  at  point  38  j  they  increase  in  fre¬ 
quency  and  move  from  the  stable  left  half  plane  across  the  imaginary  axis 
of  neutral  oscillations  (point  33)  into  the  region  of  unstable  oscillations 
in  the  right  half  plane  of  the  root  locus  diagram. 

It  is  interesting  to  compare  the  appropriate  root  locus  plot  to 
the  equivalent  locations  on  the  stability  diagram  (Fig.  13B).  Notice  in 
the  root  locus  diagram  (Fig.  13-I-b)  the  divergent  pole  at  point  37  becomes 
a  zero  root  at  point  34  which  means  that  the  E  coefficient  must  equal  zero 
at  that  point. 

On  the  stability  diagram. point  34  is  located  by  the  intersection 

Q 

of  the  value  of  =  0.0002  and  the  E  =  0  line.  By  closely  inspecting 

the  stability  diagram,  it  is  evident  that  further  increases  in  ^m  result 

M 

in  moving  upward  across  the  E  =  0  line  (point  34)  into  a  region  of  stable 
oscillations,  between  point  34  and  point  33,  on  the  stability  diagram. 

c 

Additional  increase  of  results  in  moving  across  the  R  =  0 

line  (point  33)  into  the  region  of  unstable  oscillations. 

On  the  root  locus  diagram  this  means  that  the  pole  at  point  37 

becomes  a  convergence  (point  34)  before  the  stable  oscillation  (point  38) 

G 

becomes  unstable  at  point  33.  Further  increases  in  m  on  the  root 

M 

locus  diagram  cause  the  unstable  oscillation  to  become  more  severe. 
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To  cover  the  remaining  c.g.  locations  which  are  raised  sufficiently 


to  make  m  negative  (Fig  13B),  the  root  locus  diagram  Fig._13-IV-b 

M  Q  Q 

may  be  investigated.  The  value  of  m  is  still  0.0005,  but  m  varies 

ot 

from  zero  to  negative  infinity  on  this  diagram.  The  poles  located  for 

c 

m  equal  to  zero  represent  a  convergence  (point  39),  a  divergence  (point 
M  C 

45)  and  a  stable  oscillation  (point  40  and  point  41).  As  the  gain  of  m 


is  increased  negatively,  the  poles  at  point  4  0  and  point  41  start  decreasing 
frequency  as  they  approach  the  negative  real  axis  along  the  direction  of 


the  arrows  Upon  reaching  the  entry. point  (point  47)  the  complex  pair 


becomes  two  equal,  real,  negative  roots.  As  m  continues  to  be  increased 
negatively,  the  two  equal  roots  at  point  47  split  and  move  along  the  negative 


real  axis  in  opposite  directions  with  one  terminating  at  the  zero  (point  43) 


when  is  equal  to  negative  infinity.  Meanwhile,  the  other  pole  moves 
from  47  to  the  left  on  the  real  axis  where  it  joins  with  the  original  con¬ 


vergent  pole  (formerly  at  point  39)  at  the  root  locus  break  point  (point  42}. 


At  point  42  they  become  real(equal,negative  roots  and  as  continues 


to  increase  negatively  they  break  away  from  the  real  axis  at  point  42 
to  become  a  complex  pair,  indicating  an  oscillation  of  increasing 
frequency  (point  44) . 


In  the  meantime,  there  is  always  a  divergent  pole  (point  45) 


which  becomes  increasingly  unstable  (point  46)  as  m^is  increased 
negatively. 


Again,  reference  is  made  to  Table  II  for  the  roots  of  the  quartic  for 
the  short  chord  duct  at/i=0.05,V  =  18  mph  and  the  c.g.  located  In  the 
original  position.  The  inodes  of  motion  are  two  convergences  and  an 
unstable  oscillation  with  approximately  a  ten  second  period  and  a  seven 


second  time  to  double  amplitude.  This  open-loop  response, though  unstable, 


seems  to  be  more  reasonable  than  the  others.  In  fact,  if  m  is  decreased 
by  raising  the  c.g.  ,  the  configuration  is  located  in  the  small  stable  region 


shown  on  the  stability  diagram  between  point  33  and  point  34  in  Fig.  13B. 
Actually,  this  stable  region  is  so  small  that  it  is  probably  insignificant. 


I 
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7.  Short  Chord  Duct  /i  =  0.10  V  =  36  mph 

Referring  to  the  boundaries  of  the  stability  quartic  (Fig.  14Aj; 

c 

m  changes  approximately  0.003  per  inch  of  vertical  c.g.  movement 
M  Q 

while  rn  changes  0.0001  per  inch  of  movement.  All  the  center  of 
a 

gravity  locations  appear  in  the  divergent  region  where  the  coefficients  C, 

D  and  E  are  negative,  except  for  the  lowest  location  (point  54)  where  only 
C  and  E  are  negative.  The  highest  c.g.  location  has  such  a  large  negative 

Q 

value  of  m  that  it  is  not  on  the  figure. 

M 

Starting  in  the  region  labeled 'Unstable  oscillationd1  and  decreasing 

Q 

m  along  the  line  of  centers  of  gravity,  the  E  coefficient  changes  sign 
M  Q 

at  point  52.  Slightly  below  point  52  the  axis  is  crossed  and  m  becomes 

C  ^ 

negative.  By  increasing  m  negatively  the  C  coefficient  changes 

M  Q  Q 

sign  at  point  53  and  the  location  of  m  and  m  .for  the  c.g.  located  six 

/i  a 1 

inches  below  the  originals  finally  reached  (point  54).  Raising  the  c.g. 

Q 

still  further,  the  D  coefficient  changes  sign  at  point  53  and  continues 
to  increase  negatively  as  the  higher  c.g.  locations  are  approached 
(point  56). 

Referring  to  the  stability  diagram  (Fig.  14B);  all  the  c.g.  locations 
fall  in  the  region  of  dynamic  instability.  The  points  are  below,or  on  the 
negative  side  of, the  E  =  0  line. 

The  characteristic  modes  of  motion  for  these  center  of  gravity 
locations  can  be  obtained  from  the  root,  locus  diagram. 

c  c 

Varying  negatively  at  m^  equal  to  0.001  approximately 
covers  the  region  in  which  the  centers  of  gravity  are  located.  In  Fig. 

Q 

14-IV-c  when  equals  zero,  one  pole  is  located  slightly  to  the  right 
of  the  origin  and  represents  a  divergence  (point  60).  Another  pole  is  located 
in  the  left  half  plane  and  represents  a  convergence  (point  61).  The  other 
two  poles  (point  62)  are  a  complex  pair  located  in  the  right  half  plane  and 


o 
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represent  an  unstable  oscillation.  As  the  gain  of  is  increased 
negatively,  the  pole  at  point  61  moves  to  the  right,  and  has  a  limiting 

Q 

value  at  the  zero  (point  63)  when  is  infinite.  This  pole  (point  61) 
remains  in  the  stable  region  or  left  half  plane  and  always  represents  a 


convergence  which  takes  place  with  decreasing  rapidity  as 
increased  negatively 


m  is 


The  two  poles  at  point  62  move  along  the  direction  of  the  arrows 

Q 

(point  65)  as  m  is  increased  negatively.  The  lightly  damped  unstable 
M 

oscillation  (point  62)  increases  frequency  and  becomes  a  neutral  oscillation 
(point  64)  and  then  an  unstable  oscillation  with  ever-increasing  frequency 
(point  65). 

The  predominating  mode  of  motion  for  this  quadrant  is  the  divergent 
pole  (point  60)  located  in  the  right  half  plane,  or  unstable  region,  when 

Q 

m  equals  zero.  This  pole  moves  to  the  right  (point  66)  and  causes 


the  transient  motion  of  the  machine  to  become  more  unstable  as 
is  increased  negatively. 


m 


If  the  center  of  gravity  is  lowered  considerably  more  than  six  inches, 

Q 

through  some  design  change,  m  becomes  positive.  Investigating  the 

M  Q  Q 

root  locus  using  the  extrapolated  value  of  m  »  0.001  and  m  varying 

a  /i 

from  zero  to  positive  infinity,  the  root  locus  appears  as  Fig.  14-I-c. 

C 

The  poles  for  m^ equal  to  zero,  indicate  a  convergent  mode  at  point  67, 
an  unstable  oscillation  at  point  69,  and  the  fourth  pole  (point  68)  close 
to  the  origin  is  actually  in  the  right  half  plane  (shown  by  close  analysis 

Q 

of  Fig.  14-I-h  or  Fig.  14C-b  for  m  increased  positively  to  0.001 
C  a 

while  m  equals  zero)  and  represents  a  divergence.  As  the  gain  of 

C  ^ 

m  is  increased  positively  from  zero,  the  pole  representing  a  divergence 

r 

(point  68)  moves  to  the  left  toward  the  origin  and  into  the  left  half  plane 
to  become  a  convergence  (point  71).  The  pole  at  point  67  always  moves 
to  the  left  (point  70)  and  remains  in  the  left  half  plane  as  a  convergence. 
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The  complex  pair  (point  69)  remains  unstable  and  represents  an  oscil¬ 
latory  mode  of  increasing  frequency  and  instabilityjpoint  72). 

C  C 

For  the  case  of  the  actual  c.g.  locations,  m  positive  and  m 

a  /x 

negative,  as  shown  on  the  stability  diagram,  the  open  loop  response  is 
a  serious  instability. 

C  C 

Taking  the  specific  case  of  the  values  of  m  and  m  at  the 

C  a 

original  c.g.  location  (where  m  is  admittedly  very  unstable),  the 
roots  of  the  characteristic  equation  (Table  II)  for  /i  =  0.10  V  =  36  mph 
clearly  define  the  modes  of  motion.  The  convergence  and  stable  oscilr 
lation  are  strongly  masked  by  an  instability  in  the  form  of  a  divergence 
that  doubles  in  amplitude  in  less  than  a  second  (  T  -  0.8068  seconds). 


8.  Short  Chord  Duct  Artificial  Stabilization  Added  /x  =  0.10  V  =  36 

For  the  short  chord  duct  at  fi  =  0.10,  there  is  a  dynamic  instability 
in  the  form  of  a  divergence. 

Analyzing  the  limit  cases  of  pure  rate  feedback  (T  =  0)  and  pure 
attitude  feedback  (  T  =  oo) ,  the  root  loci  appear  as  in  Fig.  15.  The  pole 
at  point  120  represents  the  instability  due  to  the  divergence  and  the 
zeros  are  located  for  increasing  K  from  zero  to  infinity.  When  T  is  equal 
to  infinity,  the  pole  and  the  zero  located  at  the  origin  cancel  each  other. 
When  T  is  equal  to  zero,  the  pole  is  located  at  negative  infinity  and  the 
zero  remains  fixed  at  the  origin. 

In  the  case  of  pure  attitude  feedback  (Fig.  15a),  where  the  pole 
cancels  the  zero  at  the  origin  (T  equal  to  infinity)  the  remaining  poles 
represent  a  stable  oscillation  (  point  115),  a  convergence  (point  117), 
and  the  divergence  (point  120).  As  K  is  increased  the  stable  oscillation 
(point  115)  continues,  but  increases  in  frequency  (point  116).  The  poles 
at  point  117  and  point  120  move  toward  each  other  until  they  become 
equal  roots  at  the  root  locus  break  point  (point  119)  in  the  right  half 
plane.  As  K  is  increased  further,  the  poles  break  away  from  the  positive 
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real  axis  (point  119)  and  become  an  unstable  oscillation,  a  neutral 
oscillation  (point  121),  and  finally  a  stable  oscillation  with  a  iim.ting 
value  at  the  zeros  (point  118)  in  the  left  half  plane. 

For  the  case  of  pure  rate  feedback  (Fig.  15b),  in  which  a  zero 
(point  136)  is  located  at  the  origin  (T  equal  to  zero),  the  same  modes  of 
motion  prevail  as  before  when  K  is  equal  to  zero.  However,  in  this  case, 
as  K  is  increased,  the  poles  at  point  130  move  toward  the  break  point 
(point  131)  and  split  so  that  one  moves  toward  point  132  and  the  other 
toward  point  133.  The  pole  reaching  point  133  joins  the  pole  formerly  at 
point  135  to  break  away  from  the  real  axis  into  the  stable  oscillation 
which  reaches  a  limiting  value  at  the  zeros  (point  134).  Meanwhile,  the 
pole  (point  137)  representing  the  divergence  can  only  become  a  zero  root 
at  point  136  when  K  is  increased  to  infinity. 

In  order,  to  study  the  effects  of  various  degrees  of  lag,  an  extra 
pole  is  situated  in  the  left  half  piano,  very  clo.^e  to  the  zero  (point  136). 
As  K  is  increased,  it  will  join  with  the  pole  at  point  135  and  break  away  as 
a  stable  oscillation  terminating  at  the  zeros  (point  134).  The  poles  at 
point  130  remain  as  stable  oscillations  which  move  awayfrom  the  real 
axis  and  continue  as  stable  oscillations  of  ever-increasing  frequency. 

As  this  extra  pole  is  moved  to  the  left,  away  from  the  zero  (point 
136),  the  former  description  of  the  loci  remains  much  the  same  except  as 
it  moves  past  the  vicinity  of  point  131  toward  point  132.  When  the 
extra  pole  moves  in  this  direction  it  pulls  the  ioci  of  point  130  closer  and 
closer  to  the  real  axis  until,  at  some  large  negative  location,  the  poles 
(point  130)  join  on  the  real  axis  similarly  to  that  shown  at  point  131 . 

The  pole  located  at  point  137  always  moves  to  the  zero  (point  136) 
and  the  basic  path  is  not  influenced  by  the  extra  pole  moving  in  the  left 
half  plane. 
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From  this  analysis  it  is  possibl  j  to  note  that  pure  attitude  feedback 
can  stabilize  this  system  {point  121  and  point  118)  for  some  value  of  gain, 
K,  whereas  pure  rate  feedback  can  only  yield  a  zero  root  (point  136). 

The  gyro-bar  with  either  aerodynamic  or  mechanical  damping 
represented  by  the  transfer  function  for  a  lagged  rate  feedback, 
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1  +  Td  ' 

should  not  be  referred  to  as  a  method  of  providing  rate  and  attitude  feedback. 
These  devices  do  not  sense  attitude  and  therefore  provide  only  rate  feedback 
with  a  phase  lag.  While  it  may  be  true  that  a  lagged  rate  feedback  provides  a 
frequency  and  phase  characteristic  similar  to  that  of  combined  rate  and 
attitude  feedback,  this  is  only  true  within  a  certain  frequency  range  so 
that,  strictly  speaking,  there  really  is  no  direct  equivalence. 

For  this  reason  the  gyro-bar  devices  may  be  able  to  diminish  the 
degree  of  dynamic  instability  of  an  aircraft  in  forward  flight  associated  with 
a  divergence  due  to  a  static  instability,  (E  <  0) ,  but  they  can  never  completely 
stabilize  the  system.  This  is  exactly  the  case  for  the  short  chord  duct  at 
jx  equal  to  0.10  and  the  e  g.  in  the  original  location.  Even  at  infinite  gain, 

K,  the  divergent  pole  (point  137)  can  only  be  made  a  zero  root  (point  136). 

An  understanding  of  the  cnaracteristic  modes  of  motion  for  configura- 

Q 

tions  of  the  Flying  Platform  with  other  combinations  of  values  for  m  and 

C  ^ 

m^  not  discussed  is  possible  by  reference  to  the  other  root  loci  provided 

for  all  the  quadrants  of  each  stability  diagram. 
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Generally,  the  unstabili  z«d,  seven-foot  diameter,  ducted  rotor 
F lying  Platform  exhibits  a  degree  of  control-fixed  dynamic  instability 
with  either  the  long  or  short  chord  duct  for  all  vertical  center  of  gravity 
locations . 

In  hovering,  both  duct  configurations  are  dynamically  unstable 

c  <*> 

due  to  the  large  positive  values  of  in  .  instability  is  characterized 

by  an  unstable  oscillation  for  all  center  of  gravity  locations  considered. 

In  forward  flight,  both  duct  configurations  exhibit  dynamic 
instability  characterized  by  unstable  oscillations  for  the  lower  c.g. 
locations  and  divergences  for  the  higher  c.g.  locations. 

For  any  of  the  configurations,  artificial  stabilization,  using 
a  rate  feedback  with  a  phase  lag  as  provided  by  a  gyro-oar,  is  able  to 
correct  a  dynamic  instability  when  it  is  of  the  form  of  an  unstable 
oscillation.  When  the  instab  iity  is  characterized  by  a  divergence  the 
lagged  rate  feedback  device  cannot  completely  overcome  the  instability. 
However,  a  combination  of  rate  and  attitude  or  pure  attitude  feedback 
can  stabilize  this  divergence. 

The  commonly  used  gyro-bars  with  aerodynamic  or  mechanical 
damping  can  only  provide  ratp  fpedhack  with  a  phase  lag  and  never  provide 
any  of  the  attitude  feedback  required  to  completely  correct  the  instability 
associated  with  a  positive  real  root  of  a  characteristic  equation.  However, 
lagged  rate  feedback  may  certainly  help  to  minimize  this  type  of  an 
instability  even  though  it  can  never  completely  overcome  it. 

The  large  tilt  angles  oredicted  for  the  faster  trimmed  flight  speeds 
will  become  a  problem  for  either  duct  configuration.  For  the  ducted  rotors 
studied,  these  excessive  tilt  angles  cannot  be  attained  because  the  pilot 
has  insufficient  control  moment  input  to  trim  the  large  moments  generated 
by  the  aircraft  for  the  higher  speeds. 

In  general,  for  given  c.g.  locations  the  moments  produced  by  the 
short  chord  duct  were  less  than  the  long  chord  duct,  but  the  moment  could 
be  reduced  for  eithr.  duct  by  raising  the  center  of  gravity. 
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Any  attempt  to  increase  the  attainable  forward  speed  by  compromising 
between  these  duct  designs  or  center  of  gravity  locations  does  not 
produce  significant  or  desirable  results.  Reducing  the  moment  by  raising  the 

c.g.  causes  the  divergent  mode  to  appear  (point  where  the  duct  becomes 

aCM 

statically  unstable  or^ —  =  E  =  0)  at  successively  lower  flight  speeds. 
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CHARACTERISTIC  EQUATIONS  OF  THE  LONG  C  fORD  DUCT 
Center  of  gravity  located,  in  the  original  position  (0.42c) 

T  -  Time  to  double  or  half  amplitude  P  -  Period  of  oscillation 


Stability  Cubic  ^overing  /l  =  V  =  0  r  =  0.3 864 

A.  +  0.31272  A  +  0.0461*25*  +  O.37872  *  0 


Roots 

A  =  -  0.81973 

T  =  O.32666  seconds 

Roots 

a2  -  0.2535  +  0.630691 

x  =  0.2535  -  0.630691 

3 

T  =  I.O563  seconds. 

P  =  3*8495  seconds 

Stability  Quartic  /i  =  0.05  V  =  l6  mph  T  =  0.3944 

A4  +  0 . 53207 x3  -  0.00007544  A2  +  0.055276 x  +  0.00868  =  0 

Roots 

A  =  -  0.13976 

T  =  1.956  seconds 

Roots 

*2  =  -  0*63364 

T  =  0.43135  seconds 

Roots 

*3  =  0.12066  +  0.288871 
x4  =  0.12066  -  0.288871 

T  =  2.265  seconds 

P  =  8.5786  seconds 

Stability  Quartic  M  =  0.10  V  =  34  mph  r  =  O.3809 

x4  +  0.62572 x3  -  0.0031x2  +  0.0282 x  +  0.010023  =  0 

Roots 

*fl=  -  0.21140 

T  -  1.2 1*86  seconds 

. 

Roots 

^2  =  -  O.66036 

T  =  0.39973  seconds 

1  Roots  , 

*3  =  0.12302  +  0.238051 

a4=  0.12302  -  0.238051 

T  =  2.1457  seconds  i 

P  =10.054  seconds 

TABLE  II 


CHARACTERISTIC  EQUATIONS  OF  THE  SHORT  CHORD  DUCT 
Center  of  gravity  located  in  the  original  position  (0.42c) 

Time  to  double  or  half  amplitude  P  -  Period  of  oscillation 


Stability  Cubic  Hovering  ^  =  V  =  0  r  =  0. 33926 

x3  +  0. 34493  x2  +  0.014260 x  +  0.073105  =  0 

Roots 

X  =  -0-55575 

T  =  0.42304  seconds 

Roots 

X2  *  0.1054  +  0. 34684i 
*  A-3  =  0.1054  -  0. 346841 

T  =  2.2306  seconds 

P  =  6.1459  seconds 

Stability  Quartic  ^  =0.05 

V  =  18  mph  r  =  O.35I 

X4+  0. 45072 x3  +  0.04o44  x2 

+  0.02577 x  +  0.0013625  *  0 

Roots 

- 

=  -0.05511 

T  =  4.4l4  seconds  - 

Roots 

X2  -0.46847 

T  =  0.5192  seconds 

Roots 

x3=  +  0.03643  +  0.226821 

T  =  6.677  seconds 

+  0.03643  -  0.226821 

P  =  9-7231  seconds 

TABLE  III 
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OVERALL  VALUE  Of  LIABILITY  DERIVATIVES 
FOR  VARIOUS  C.  G.  LOCATIONS 
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FIGURE  1  FLYING  PLATFORM 


FIGURE  2  FULL  SIZE  DUCT  PROFILES  FOR  WIND  TUNNEL  MODEL  TESTS 
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All  notation  shown  is  positive  except  flight  path 
angle  to  horizon  ,  7  . 


Velocity  Ratio 


FIGURE  5  TYPICAL  TIME  RESPONSE  TO  A  SMALL  INPUT  IN  HOVERING 


Velocity  Rati 


FIGURE  6  TYPICAL  TIME  RESPONSE  TO  A  SMALL  INPUT  IN  FORWARD  FLIGHT 


FIGURE  7  ROOT  LOCUS  S.C.D.  AND  L.C.D.  HOVER 
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FIGURE  8  STABILITY-  DIAGRAM  FOR  HOVER  C  >  0 
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FIGURE  9  ROOT  LOCUS  S.C.D.  HOVER-ARTIFICIAL  STABILIZATION  ADDED 
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BOUNDARIES  OF  THE  STABILITY  QUARTIC  S.C.  D.  =  0.05  V  =  18MPH 
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EVALUATION  OF  STABILITY  DERIVATIVES 

The  stability  derivatives  were  evaluated  from  data  obtained  by 
wind  tunnel  tests  conducted  at  the  David  Taylor  Model  Basin.  A  two  foot 
diameter,.  2/1  scale,  model  of  a  ducted -rotor  Flying  Platform  was  tested 
in  the  DTMB  low-speed  wind  tunnel  (Ref.  2).  From  the  large  variety  of 
test  runs  made,  portions  of  the  data  incorporating  different  duct  designs 
were  selected  for  the  purpose  of  obtaining  stability  derivatives.  The  duct 
designs  chosen  were:  one  of  relatively  long  chord  length  with  a  large 
leading  edge  radius,  the  other  of  shorter  chord  length  and  smaller  leading 
edge  radius  (Fig.  2).  The  latter  was  chosen  because  it  represented  a  con¬ 
siderable  design  change  from  the  large  leading  edge,  long  chord  duct. 

The  tests  were  conducted  in  the  large  settling  chamber  of  the 
wind  tunnel.  The  balance  was  designed  so  that  forces  were  measured 
perpendicular  and  parallel  to  the  rotor  shaft  and  the  moments  were 
measured  about  the  supporting  structure. 

# 

The  moment  was  arbitrarily  transferred  to  the  so-called  quarter 
chord  point  ot  the  duct.  The  forces  and  moments  so  measured  were  T, 

H  and  Mc/4  (Fig.  3). 

All  runs  were  made  at  constant  velocity  ratios  as  the  angle  of  attack 
of  the  duct  was  varied.  The  tests  were  conducted  principally  at  only 
three  different  velocity  ratios,  making  accurate  determination  of  derivatives 
with  respect  to  pi  somewhat  difficult. 

When  summing  all  the  forces  and  moments  acting  on  the  machine,  it 
was  necessary  to  include  the  drag  force  and  moment  caused  by  the  pilot 
in  addition  to  the  forces  and  moments  measured  in  the  wind  tunnel,  (Fig.  3) 
This  drag  force,  directed  opposite  and  parallel  to  the  relative  wind,  was 
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taken  as  acting  through  an  assumed  drag  center  of  the  pilot.  In  this  report, 
the  drag  center  for  the  typical  pilot  is  located  3^  feet  above  the  pilot's 
feet . 

Various  methods  are  available  to  determine  the  drag  coefficient  or 
equivalent  flat  plate  area  of  the  pilot.  A  report  published  on  the  subject  of 
airloads  on  human  beings  (Ref.  9),  determines  the  drag  coefficient  of 
varo.  s  men  in  standing  positions.  For  the  specific  case  of  a  six-foot  , 

200  pound,  clothed  man,  standing  erect,  with  arms  and  feet  together,  the 
equivalent  flat  plate  area  was  determined  as  approximately  nine  and  one 
half  square  feet.  No  attempt  was  made  to  modify  this  value  due  to  the 
pilot  stance  with  arms  extended  and  feet  apart,  or  to  compensate  for  ad¬ 
ditional  components,  such  as  the  cage  around  the  pilot,  instrumentation, 
fuel  tanks,  pilot  helmet  and  equipment.  It  was  felt  that  attempts  to  make 
accurate  modifications  on  the  drag  of  the  pilot,  in  any  way,  would  not 
be  worthwhile  unless  investigations  were  made  of  the  air  flow  to  which 
the  pilot  is  exposed.  Therefore,  for  the  purposes  of  this  report,  the 
value  of  the  equivalent  flat  plate  area  of  the  pilot  is  assumed  to  be 
nine  and  one  half  square  feet. 

The  original  center  of  gravity  position  for  the  complete  machine 
at  normal  gross  weight  was  arbitrarily  located  at  the  pilot's  stand-on 
platform  or  0.42c  above  the  top  of  the  duct.  Additional  vertical  center 
of  gravity  positions  were  selected  above  and  below  this  original  reference 
point  (Fig.  16).  All  moments  were  resolved  about  these  center  of  gravity 

locations  and  the  stability  derivatives  evaluated  for  each  case.  (Table  III). 

* 

The  important  requirement  in  the  evaluation  of  the  stability  deri¬ 
vatives  is  that  they  be  calculated  at  each  initial,  steady-state,  equili¬ 
brium  trim  condition  of  the  aircraft.  For  this  steady  flight  condition,  the 
sum  of  all  the  aerodynamic  forces  in  the  drag  direction  must  be  equal  to 
zero  and  the  sum  of  the  forces  in  the  lift  direction  must  equal  to  the  gross 
weight  of  the  aircraft. 
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Fig.  20  Equilibrium  Diagram 


In  Fig.  20  where  or  is  a  small  angle,  the  summation  of  forces 
in  the  drag  direction  is; 


IE 
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=  -  Cr  sm  otf  -  cos  «r  -  C//f  *=  O 


This  equation  may  be  solved  to  determine  the  variation  of  the 
equilibrium  angles  of  attack  with  forward  speed. 

Actually,  this  was  accomplished  graphically  by  plotting  the 
horizontal  force  versus  angle  of  attack  for  the  various  forward  speeds 
(Fig.  21). 
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Values  of  equilibrium  angles  plotted  versus  speed  are  shown  for 
both  ducts  in  Fig.  18  on  page  117. 

For  vertical  equilibrium  in  the  initial  steady-state  level  flight 
condition,  (when  the  sum  of  the  forces  in  the  lift  direction  equals  the 
gross  weight)  it  is  necessary  to  determine  the  variation  of  rotor  angular 
velocity  with  forward  speed  from: 

L  =  w=  (+ Cr  **  *  -  e„  s*c<r)  j>  7rR2(nR)z  <64) 

where  R  =  3.5  feet,  W*  550  pounds  and  CTT,  and  a  are  obtained 

H  T 

from  the  wind  tunnel  data  for  the  particular  equilibrium  angle  and  velocity 
ratio . 

Then,  for  vertical  equilibrium  the  actual  flight  speed  may  be 
determined  from: 

a  ClR 

'  (65) 

Actual  values  of  &R  at  the  various  velocity  or  tip  speed  ratios 
are  listed  for  both  ducts  in  Table  IV. 
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TABLE  IV 

ROTOR  TIP  SPEED  (it/sec)  versus  yi 


,  Velocity  Ratio 

o 

ii 

St 

M  =  0. 05 

M  3  o.io 

1  &R  Short  Chord  Duct 

555 

537 

525 

OR  Long  Chord  Duct 

- - - L 

4Z3 

478 

495 

When  evaluating  the  stability  derivatives,  the  forces  and  moments 
were  put  into  coefficient  form  and  plotted  versus  the  angle  of  attack 
(measured  in  radians)  and  ji  (dimensionless  velocity  ratio). 

Typical  sketches  in  the  evaluation  of  the  derivatives  are  shown 
in  the  following  Figures. 


The  short  vertical  lines,  on  the  constant  ^  curves ,  indicate 
the  previously  calculated  equilibrium  angles  for  each  duct  at  each 


*  0.  •  upv  i 


The  stability  derivative  CT '  Is  obtained  by  taking  the  slope  at 
the  equilibrium  angle  for  the  particular  ft. 

Occasionally  it  was  necessary  to  extrapolate  in  order  to  obtain 
the  hovering  values  (jt  •  0)  of  derivatives  taken  with  respect  to  ft. 

This  can  be  avoided  if  more  values  at  zero  alpha  are  obtained  at  the 
higher  ft  or  by  additional  data  at  lower  velocity  ratios. 

Values  of  T  obtained  from  the  data  are  listed  in  Table  VI, 

TABLE  VI 


T  versus  u 

M  * 


w 

Velocity  Ratio 

M  =  0. 05 

/2  =  0.10 

C 

T  Short  Chord  Duct 

0 

0.044 

0.  016 

c 

T  Long  Chord  Duct 

0 

-0.004 

-0. 044 

For  CH  : 


a 


Equilibrium  angles  of  attack  for  /x  =  0.15  were  usually  so 
large  that  they  were  not  shown  on  the  figures. 


Values  of  H  are  listed  in  Table  VII 
a 


TABLE  VII 
C  versus  /i 


Velocity  Ratio 


C  Short  Chord  Duct 

ri 

a 


Long  Chord  Duct 


^=0 

IL  =  0.05 

0 

0 

0 

0.0005 

0.0015 


m 


— T*  ' 


Cress  plotting  for  C  : 


a  =  0,  -0.1,  -0.2 


A  =  equilibtium  a 


Fig.  25  Sketch  of  C..  versus  n 

H 


Values  of  G  are  listed  in  Table  VIII. 
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TABLE  V 


C  versus 
H 


Velocity  Ratio 

n=o 

/i  =  0.05 

/i  =  0. 10 

c 

H  Short  Chord  Duct 

/* 

0.044 

0.078 

0. 108 

c 

H  Long  Chord  Duct 

/i 

0,118 

0. 118 

0.136 

equilibrium  a 


Velocity  Ratio 


m  Short  Chord  Duct  I  0 
a 


m  Long  Chord  Duct 
a 


0.05 

0. 10 

0.00035 

0.0023 

0.0031 

0.0046 

Velocity  Ratio 


m  Short  Chord  Duct 


0. 055(ext. ) 


m  Long  Chord  Duct 
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Values  for  H 


ft,  z 


H/R 


d  (^0) 
,  and  r  are  listed 


using  the  appropriate  values 
in  Table  XI. 


of 


TABLE  XI 


C  versus  y 

Hd  (AS) 


Velocity  Ratio 

1!  > 

o 

y  =  0.05 

y  =  0.10 

CH  , ,  ..Short  Chord 
d(a9)  Duct 

0.000275 

0.000486 

0.000673 

c 

H , ,  .  Long  Chord 

d  (  A  0)  _ 

v  '  Duct 

0.000735 

0.000735 

0.000848 

In  the  equations  of  motion  the  drag  terms  (Terms  5,  6  and  7 
on  page  19)  are  developed  as  a  result  of  the  drag  of  the  pilot  and  may  be 
evaluated  directly.  Their  form  in  the  equations  may  also  be  seen  in 
Fig.  20  where  the  coefficient  form  of  this  force  (acting  parallel  to  the 
relative  wind)  is  y  C^'  /2  (Eq.  13). 

Using  the  small  angle  assumptions,  the  pilot  drag  force, 

2 

ii  CnK/2  is  resolved  into  components  along  the  body  axes.  The  derivatives 
D 

are  taken  with  respect  to  yand  a  and  must  be  evaluated  at  the  initial 

condition,  a  -  u.  and  a  =  a  . 

2  ° 

X  component:  y  C^1 
2 

2  ^  , 

Z  component:  y  a 


Taking  the  derivatives: 


X  axis: 

z  dn 

CD  a 


(*n  =  jicD' 
(Jl2)  =  o 


o 


-  “oV 


are  obtained  by  multiplying  these 


terms  by  their  respective  perturbations 


H-force  contribution 


Thrust  force  contributions 


APPENDIX  II 
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The  remaining  stability  derivative  to  be  evaluated  is  the  over-all 

c 

pitch  damping  derivative  mrf  (Eq.  .21g), 


c  =  c 

d(^6)  d(A  Q) 


-  2 


/?  * 


(21g) 


The  total  damping  moment  in  pitch  for  the  unstabilized  Flying 
Platform  includes  contributions  from  three  sources,  namely: 

1.  The  pitch  damping  moment  due  to  the  change  in  H-force  with 


pitch  rate,  hL.C 


H 


.  a.(±e) 


R  XAd  {  a  Q) 

2,  The  pitch  damping  moment  of  the  rotor  due  to  a  pitch  rate. 

3.  The  pitch  damping  moment  of  the  duct  due  to  a  pitch  ■■rate. 

Contributions  (2)  and  (3)  taken  together  are  equai  to  the  first  term 
C 

m  . , .  in  Eq.  21g. 


d(A<?) 


Contribution  (1)  is  obtained  by  solving  Eq.  68  for 


d(  A  6) 


L 


Hd(*e) 


-  C,  .  I±L  .  1 

r  siy 


The  values  of  H 


12  ,  Rand  t  are  all  known.  Then 


the  pitch  damping  moment  contribution  due  to  change  of  H-force  with 
pitch  rate  is: 


(68) 


?Ltt  .  CV,  >cL(&Q) 
R  d(*6) 


(Term  17  page  21) 


Contribution  (2)  is  obtained  by  using  the  expression  developed 
in  the  Appendix  of  Ref.  4  for  the  moment  generated  by  the  rotor  due  to 
a  pitch  rate. 


pitch  rate  about 
Y  axis  at  top  of 
duct. 

duct  exit  area 
duct  chord  length 
slipstream  velocity 
at  duct  exit. 


us  of  gyration  for  the  duct  exit  area 


The  duct  exit  velocity  is  equal  to  (See  Appendix  lit) 


Therefore  the  pitch  damping  moment  in  coefficient  form  is 


(*e)j>7rtfgR\ cy1  r  p 
"pvR'tflRTR  L  + 


The  pitch  damping  derivative  i 


d(*e)  *  * 

where  w  is  equal  to  the  aircraft  relative  densit 
For  the  long  chord  duct  when  R  =  3.5  ft 
The  expression  for  contribution  (3)  is; 


over 


The  sum  of  the  three  moment  contributions  gives  the  total  theoretical 
pitch  damping  moment  used  in  this  Report. 

It  is  recognized  that  this  theoretical  pitch  damping  moment  is  probably 
too  large,  since  it  is  not  very  likely  that  all  three  contributions  provide  their 
full  theoretical  value  for  the  cases  of  the  relatively  short  chord  ^ucted  rotors 
studied.  In  any  event,  this  total  value  appears  to  be  small  compared  to  the 
amount  required  tc  produce  desirable  stability  characteristics.  Therefore, 
the  stability  characteristics  together  with  the  effects  produced  by  additional 
damping  were  investigated  and  are  liscussed  in  this  Report. 
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APPENDIX  III 


The  expression  for  duct  exit  velocity  may  be  developed  from 

the  so-called  simple  momentum  theory  of  ducted -propellers  by  utilizing 

the  common  assumption  that  the  velocity  of  the  air  at  the  duct  exit 

is  equal  to  the  final  slipstream  velocity  (V  .  =  V  ). 

exit  oo 

The  analysis  is  restricted  to  the  inviscid,  incompressible, 
static  case  and  utilizes  the  ideal  ducted  propeller  combination  with  an 
“actuator  c'isk"  and  a  simple  cylindrical  shroud. 


The  velocity  V  is  the  same  on  both  sides  of  the  actuator  disk 
because  of  the  requirement  of  continuity. 

The  total  thriist  is: 


Utilizing  the  assumption  of  no  slipstream  contraction  for  a 
ducted  pi  speller  Ap  =  A^,  *  A  . 

By  continuity: 
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